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ABSTRACT 


This report describes the results of a 2 -year study of the 
mechanism of spall failure in the HPOTP bearings. The objective 
was build a foundation for detailed analyses of the contact life 
in terms of; (i) cyclic plasticity, (ii) contact mechanics, (iii) 
spall nucleation and (iv) spall growth. Since the laboratory 
rolling contact testing is carried out in the 3-ball-rod contact 
fatigue testing machine, the analysiss of the contacts and 
contact lives produced in this machine received attention. 

The analysis of previous cyclic stress-strain hysteresis 
loop measurements of 44 0C steel was refined to account for the 
plasticity of the fillet regions. In addition, the hysteresis 
loop shapess of the hardened 7075 aluminum alloy were measured. 
In both cases the elastic-linear-kinematic-hardening-plastic 
(ELKP) loop parameters were evaluated. Elasto-plastic, finite 
element analyses of the repeated, 3-dimensional, frictionless, 
rolling contact produced in the 3-ball-rod testing machine at 
Hertzian pressures of p G = 2.4, 4.0 and 5.4 GPa (for 440C steel) 
and 1.25 GPa (for hardened aluminum) were carried out using the 
appropriate ELKP-loop parameters. These calculations were also 
extended for 440C steel properties to 3-dimensional rolling-plus- 
sliding and to the 2-dimensional (line contact) thermal- 
mechanical coupled rolling-plus-sliding with frictional heating. 
The results of calculations are compared with observations of 
aluminum rods subjected to contact under these conditions. 

Rolling contact tests of the 440C steel and the hardened 
aluminum were performed in the 3-ball-rod testing machine with 
smooth and roughened balls. Efforts were made to evaluate the 
effects of retained austenite. A series of tests were performed 
on the 440C samples with small « 100 /xm indentations in the 
running track which make it possible to locate and follow the 
progress of spall nucleation and growth. The results define the 
spall nucleation- and spall growth-component of the contact life 
of the 440C steel over a range of the contact pressures. They 
also provide evidence for a threshold pressure for crack growth. 
In addition, the 3-dimensional features of the spall were studied 
by a novel replicating technique and metallographic sections of 
the spalls. 

The contributions to the fracture mechanics crack growth 
driving force for surface breaking cracks arising from the 
Hertzian stresses, surface irregularities, fluid in the crack 
cavity, centrifugal stresses and thermal stresses are reviewed 
and the results of 2- and 3-dimensional analyses are compiled and 
compared. New calculations for the Bower model of a 2-dimensional 
surface breaking crack with fluid pressure in the crack cavity 
are pressented. A numerical expression for the Mode I crack 
driving force for a 3-dimensional crack with fluid in the crack 
cavity is used to calculate the spall growth component of the 
contact life. These calculations are compared with and are in 
reasonable agreement with the measurements of spall growth. 
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1. INTRODUCTION 
1 . 1 Background 

Spall failures limit the life and performance of the HPOTP 
bearings in the space shuttle main engines. While the general 
features of spalling are known, detailed quantitative analyses of 
the contact life in terms of contact geometry, loading and basic 
material properties have not been developed. As a consequence, 
efforts to improve life proceed by trial and error and are close- 
ly tied to uncertain and time consuming laboratory testing — 
uncertain because the laboratory tests do not reproduce the 
service conditions. 

The analysis of contact life is complicated by the existence 
of 2 distinct failure modes that result from 2 sources of contact 
plasticity (Hahn and Rubin, 1990) : 

■ Subsurface Originated Spall Failure . This is caused by 
the translating, Hertzian pressure pulse which produces peak 
amounts of cyclic plasticity and damage in an annular layer at a 
depth, z » 0.5 w (2w is the contact width) below the surface. 

■ Near-Surface Originated Spall Failure . This is produced 
by the cyclic plasticity and damage produced just below the 
running track at depths, z ~ 1 /xm to 50 jum, caused by several 
sources: (i) stationary pressure spikes produced by surface 
irregularities, asperities, debris dents, etc, ( i i ) tractions 
arising from the sliding of the contact with friction, and (iii) 
thermal stresses from frictional heating. 

The distinction is important because both the nature of the 
cyclic deformation and damage leading to spall nucleation, and 
the mechanism of growth of the spall are different for the 2 
modes. Both the observations of spalls in the HPOTP bearings 
(Bhat and Dolan, 1982) and the results of this study support the 
view that the failures of the 440C steel bearings are of the 
near-surface mode. 

Work aimed at improving the performance of the HPOTP bear- 
ings can benefit from an analysis of the near-surface spall 
failure mode composed of the following elements: 

(i) Cyclic Plasticity . The definition of the continuing, 
near-surface cyclic plasticity. This is governed by the shape 
of the cyclic stress-strain hysteresis loop of the steel. While 
the 4 4 0C loop shapes were measured in a previous NASA-supported 
study (Kumar et al., 1987), subsequent work revealed that the 
analysis of the measurements requires improvement. The 
connections between the material microstructure and the stress- 
strain properties must be clarified. 
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(ii) Contact Mechanics . The definition of the cyclic 
plasticity must also draw on the mechanics of contact. This can 
be treated using elasto-plastic finite element methods. In the 
previous work (Kumar et al., 1987) the present authors have 
devised finite element analyses of 2-dimensional (line contact) 
rolling-plus-sliding, and rolling-plus-sliding with heat 
generation (Kulkarni et al., 1991). These methods need to be 
extended to the 3-dimensional contact, refined and critically 
tested. i ' 

(iii) Crack Nucleation . The rates with which the continuing 
cyclic plasticity lead to the accumulation of damage and crack 
nucleation must be formulated. The work to accomplish this is 
in its early stages (Keer et al., 1986) and can benefit from 
experimental determinations of the nucleation component of the 
contact life and the factors that influence it. 

(iv) Spall Growth . The rate of spall growth must be 
evaluated. This is governed by the fracture mechanics driving 
force and the steel's Mode I, II and III, da/dN- AK 
characteristics. The cracks produced by the near-surface mode 
become surface breaking at an early stage of their life. As a 
result, the driving force is amplified by the pressure of 
lubricant fluid forced in the crack cavity (Bower, 1988) . In the 
absence of lubricant, thermal stresses arising from frictional 
heating enhance the crack driving force (Goshima and Keer, 1990) . 
The analysis of the growth life calls for calculations of the 
driving force for small surface-breaking cracks containing fluid 
pressure as well as measurements of relevant da/dN- K 
properties of the material. Finally, there is a pressing need 
for measurements of the spall growth component of the contact 
life that can be used to test the reliability of the fracture 
mechanics method. 

1.2 Summary 

This report describes the results of a 2-year follow-on to 
an earlier NASA-supported study on the mechanism of spall 
failure in the HPOTP bearings (Kumar et al., 1987). The 
objective was to build a foundation for detailed analyses of the 
contact life along the lines described above. Since much of the 
laboratory rolling contact testing is carried out in the 3-ball- 
rod contact fatigue testing machine, the analysis of the contacts 
and contact lives produced in this machine received attention. 

The following tasks were undertaken: 

(i) Cyclic Plasticity . The analysis of the previous 
cyclic, stress-strain hysteresis loop measurements of 440C steel 
was refined to account for the contribution of plasticity in the 
fillet regions. In addition, the stress-strain hysteresis loop 
shape of the hardened 7075 Al alloy was measured and the elastic- 
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linear-kinematic-hardening plastic- (ELKP-) loop parameters were 
evaluated. These studies are described in Sections 2.3 and 2.4. 
The results are incorporated in 3 -dimens iona 1 , elastoplastic, 
finite element analyses of the 3-ball-rod testing machine contact 
conditions described in Section 3. 

(ii) Elasto-Plastic Finite Element Analyses . Elasto- 
plastic, finite element analyses of the 3-dimensional rolling 
produced in the Federal Mogul/Bowers/NTN, 3-ball-rod testing 
machine at several contact pressures were carried out using both 
the 440C steel and 7075 aluminum ELKP loop parameters (see 
Sections 3.3 and 3.4). The calculations are compared in Section 
3.4 with experimental observations on aluminum rods (tested in 
the 3-ball-rod testing machine) that characterize the size and 
shape of the 3-dimensional contact cyclic plastic zone and the 
propensity for subsurface crack nucleation and growth. In 
addition, the finite element calculations were extended to 3- 
dimensional rolling-plus-sliding (Section 3.3) and to the 2- 
dimensional (line contact) thermal-mechanical coupled problem for 
440C steel, ELKP loop properties (Section 3.5). 

(iii) Fracture Mechanics Analyses of Surface Breaking 
Cracks . Contributions to the crack growth driving force from 
the Hertzian stresses, surface irregularities, fluid in the crack 
cavity, centrifugal stresses and thermal stress are reviewed and 
the results of different 2-dimensional and 3-dimensional analyses 
compiled and compared in Section 4.4. New calculations for the 
Bower model of a 2 -dimensional surface breaking crack with fluid 
in the crack cavity were carried out and are presented in Section 
4.3. The implications of the driving force values with respect 
to the threshold crack size were examined (Section 4.5). 

(iv) Retained Austenite . Efforts were made to 
evaluate the effects of retained austenite in 440 C steel on the 
contact life and these are described in Section 5.3. 

(v) Measurements of the Nucleation and Growth Lives . 
The separate contributions of nucleation and growth components of 
the contact lives of 440C steel rods tested in the 3-ball-rod 
testing machine were measured. These experiments also reveal 
effects of surface roughness and a threshold crack size for spall 
growth. This is presented in Sections 5.4 and 5.5. 

(vi) Characterization of the 3-Dimensional Spall . 

A novel replication technique was developed to find the three 
dimensional features of the spalls, this technique is reported in 
Section 5.6. This section also reports the results of 
metallographic investigation of the main geometric features of 
the spall. 

(vii) Evaluation of the Fracture Mechanics Driving Force 
for Spall Growth . Calculations were performed to determine the 
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Mode II crack driving force under different conditions, as 
reported in Section 4.3. These results were compared with other 
works which accounted for thermal loading and three dimensional 
effects, and are reported in Section 4.4. The critical, i.e. 
threshold, crack sizes were evaluated in light of the above 
mentioned results in Section 4.5. 

(viii) Analysis of the growth life . An attempt is made to 
compare the experimentally observed lives with those predicted 
based upon the results drawn in (vii) . The conclusions are 
reported in Section 4.6. 
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2. CYCLIC STRESS-STRAIN PROPERTIES 
2 . 1 Background 

The continuing cyclic plasticity, which accompanies rolling 
contact when the Hertzian contact pressure exceeds the shakedown 
limit, damages the material and may ultimately lead to spall 
nucleation (Hahn and Rubin, 1990) . The plasticity may also 
assist the process of cyclic crack growth (Bastias, 1990) . Both 
the value of the shakedown pressure and the amounts and distribu- 
tion of cyclic plasticity that occurs above shakedown depend on 
the material's resistance to plasticity (Merwin and Johnson, 
1963, Bhargava et al., 1985, 1990, Hahn et al., 1987, Hahn and 
Rubin 1990) . For the case of repeated contacts which produce 
essentially fully reversed cyclic plasticity, the resistance is 
given by the shape of the stress-strain hysteresis loop (Hahn et 
al., 1990). The constitutive relations that describe the loop 
shape must be incorporated into the finite element models of 
contact described in Section 3. 

In the past, it has been common practice to treat the cyclic 
plasticity as isotropic and elastic-perfectly-plastic (EPP) . The 
loop produced by this highly idealized behavior, shown 
schematically in Figure 2.1a, can be described by 2 parameters: 
the elastic modulus, E, and the shear yield strength, k. Figure 
2.1b illustrates that the loop shapes of 440C steel and other 
bearing steels are not approximated by EPP-behavior . The real 
loops display rapid strain hardening and kinematic behavior 
(Hahn et al., 1990). To improve the analyses, the authors have 
devised a bilinear, 3-parameter elastic-linear-kinematic- 
hardening-plastic (ELKP) representation of the loop illustrated 
in Figure 2.2a (Hahn et al., 1987, Hahn et al. , 1990). The loop 
parameters: the elastic modulus, G, the kinematic shear yield 
strength, k^, and the plastic modulus, M s , are defined in Figure 
2.2a. The relations between these parameters and the 
conventional parameters defined in Figure 2.2a, including the 
stress amplitude, o a , and the energy dissipated (loop area), u', 
are given in Table 2.1. 

In an earlier report (Kumar et al., 1987) the authors 
described the results of cyclic torsion tests performed on 
hardened 440C steel which were used to evaluate the ELKP 
parameters. These analyses assumed that the contribution of the 
fillet region could be neglected. Subsequent work showed this to 
be a poor approximation and a method for accounting for the 
plasticity of the fillet was devised (Hahn et al., 1990). This 
method has been employed here to analyze the measurements 
reported earlier so as to provide more reliable values of the 
parameters for 440C steel. 

In addition, cyclic torsion tests were performed on hardened 
7075 aluminum, to establish the ELKP-parameters for finite 
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Figure 2.1 Shear stress-shear strain hysteresis loops: (a) 

the loop for idealized, isotropic, elastic- 
perfectly-plastic (EPP) behavior, and (b) the 
loops displayed by 440C steel after N=15 and N=250 
stress cycles. While the EPP-loop is drawn so 
that its 0. 035%-offset, shear yield strength 
corresponds with that of the N=15 loop of the 440C 
steel, it is clear that the EPP loop does not come 
close to representing the cyclic stress-strain 
behavior of the steel. 


ORIGINAL PAGE IS 
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Figure 2.2 Schematic of the bilinear, 3-parameter, elastic- 

linear-kinematic-hardening-plastic ( ELKP ) 
representation of the hysteresis loop: (a) 
conventional form employed for 440C steel and (b) 
special form employed for 7075-T6 aluminum to 
accommodate the differences in the elastic modulus 
in tension and compression. The 3 ELKP-parameters 
are: (for tension-compression) the elastic 
modulus, E, the kinematic yield strength, a^, and 
the plastic modulus (G, k^, and M s for torsion) . 
The relations between these parameters and more 
conventional parameters are given in Table 2.1. 


3! M 
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Table 2.1 Relationship Between the ELKP Loop Parameters: 

Ej Q]r . M. and the Conventional Properties of the 

Loop . 


a e = a k - [ E M / (E - M) ] £eP/2 ± a m 

a OfC (0.02%) = a k - [EM/ (E - M) ] (AeP/2 - 0.0002) ± 
a a = a k + [EM/ (E - M) ] AeP/2 
U' = 2 a k eP 


a k 

°e 

^oc ( 0 . 02-s ) 


'm 


AeP/2 


U' 


Young's Modulus 

Kinematic Hardening Modulus 

Kinematic yield strength 

Cyclic elastic limit (measured from zero stress) 
Cyclic, 0.02% offset yield strength (measured from 
zero stress) 

Stress amplitude 
Mean stress 

Plastic strain amplitude 

Per cycle plastic work or loop area 
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element analyses of rolling contact for this material for reasons 
mentioned in Sections 3.4. 

2.2 Experimental Procedures 

The Cyclic torsion tests were carried out on hollow cylin- 
drical 4 4 OC using the procedures described in a previous report 
(Kumar et al., 1987). The cyclic torsion tests were plastic 
strain amplitude controlled. The torque-rotation loop was 
recorded, converted into a shear stress-shear strain hysteresis 
loop, the ELKP-parameters : the kinematic shear yield strength, 
k k , the plastic modulus, M s , described, and the equivalent 
tensile values, a K and M, were evaluated (Hahn et al., 1990). 
Unlike the analysis employed in the previous report (Kumar et 
al., 1987) the procedure accounted for the plastic contribution 
of the filler regions, which is significant (Hahn et al, 1990). 

The studies performed on the hardened 7075-T6 1 aluminum 
alloy were carried out in uniaxial, push-pull fatigue. The test 
bar is illustrated in Figure 2.3. The average hardness of the 
as-heat treated samples was HRB-87 (HK-170) . The test pieces for 
the 3-ball-rod testing machine, referred to in Section 3.4, were 
machined from the same stock and heat treated in the same way. 
The shapes of the cyclic stress-strain hysteresis loops obtained 
under conditions of constant strain amplitude were measured at 
room temperature at a frequency of f=0.75 Hz with a servohydrau- 
lic testing machine by manually adjusting the stress amplitude. 
The axial strain was measured with an extensometer attached to 
the gage section and the hysteresis loops were periodically 
recorded and analyzed with a high speed data acquisition system 
programmed to evaluate the ELKP-parameters, and M. The 
details of the procedure are similar to the one used to analyze 
the axial torsion tests mentioned above (Hahn et al., 1990). 

In the case of the 7075-T6 aluminum, the slope of the 
elastic portions of the loop on the tensile side, E = 68.6 GPa, 
is significantly lower than the slope on the compression side of 
the loop, E = 71.5 GPa, as shown schematically in Figure 2.2b. 
These values do not change with number of stress cycles. When 
the difference in the values of the slopes was recognized and 
accounted for, as in Figure 2.2b, the values of ex* and M for the 
tensile and compressive portions of the cycle were in close 
agreement. To simplify the treatment of the constitutive 
relation, the averages of the values of E, a ^ and M obtained 
from the tensile and compression part of the cycle are quoted and 
were inserted in the finite element analyses reported in Section 


1 Composition and heat treatment of the 7075 aluminum alloy: 
Zn-5.6, Mg-2.5, Cu-1.6 and Cr-0.23; the alloy was solution 
treated at 870 °F, spray quenched in water and aged at 250°F for 
24 hours. 



Figure 2 . 3 


Axial fatigue test specimen used to measure 7075 
T6 hysteresis loops. All dimensions are i 
meters . 
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3.4. The difference between the tension and compression modulus 
observed here is similar to that quoted in the ASM Metals 
Handbook (Vol.2, 9th Edition, 1979). The origin of the 
difference is not known to the authors. It is possible that a 
large number of the small, 5 /im diameter particles, in the alloy 
may rupture during the tensile portion of the first few stress 
cycles producing small cracks that reduce the stiffness in 
tension but not in compression. 

2.3 Cyclic Stress-Strain Properties of 440C Steel 

Results of the hysteresis loop measurements are summarized 
in Figures 2.1b, 2. 4-2. 7 and in Table 2.2. Examples of the 
hysteresis loops displayed by 440C steel samples are reproduced 
in Figure 2.1b. The 0.035% shear strain offset (equivalent to an 
0.02% tensile strain offset) is marked by the short horizontal 
tick in Figure 2.1b. The application of a mean stress produces 
a not-fully reversed loop — more plasticity in one direction 
than in the other. However, the non-reversibility 2 diminishes 
with increasing number of cycles and effectively disappears after 
N'10 3 cycles as shown in Figure 2.4. This is viewed as evidence 
that the cyclic behavior closely approaches kinematic behavior 
after N > 10 3 cycles. The systematic decrease of the conven- 
tional cyclic yield strength, Oqq, with increasing plastic strain 
amplitude, and the negative value of this quantity at the largest 
strain amplitude, are also incompatible with isotropic behavior 
and consistent with kinematic behavior. 

The variation of the kinematic yield strength, a^, with 
numbers of cycles is shown in Figure 2.5. All of the samples 
display noticeable cyclic hardening — increases in with 
number of cycles — during their relatively short lives: 300 < N 
< 600, and are still hardening after N = 600 cycles. The near- 
end-of-life values 3 reported in Table 2.2 represent a lower bound 
estimate of the kinematic yield strength appropriate for large 
numbers of cycles, e.g. N = 10^ to 10^. These a^-values are 
about 15% lower than the values reported earlier which were 
obtained neglecting the plasticity in the fillet regions. The 
values of the plastic modulus, M, in Figure 2.6 show virtually no 
dependence on the number of cycles. These values are about 2% 
smaller than those reported earlier. 


2 Non-reversibility is defined as: NR = i-eP f /eP r , where 
ePf is the plastic strain increment in the forward direction and 
eP r is the increment in the reverse direction. Non- 
reversibility is associated with isotropic cyclic behavior and is 
absent for kinematic cyclic behavior. The decay of the non- 
reversibility is symptomatic of the approach to kinematic behavior. 

3 Values for the largest number of stress cycles that 
appear unaffected by the nucleation and growth of a crack. 
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Figure 2.4 Variation of the non-reversibility of the 440C 

steel hysteresis loop produced by a mean shear 
stress, r m = 200 MPa, with number of cycles. 


KINEMATIC YIELD STRESS, 



NUMBER OF CYCLES, N 


Figure 2.5 


Variation of the kinematic yield strength of 440C 
steel with number of cycles. 
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Variation of the plastic modulus of 440C steel 
with number of cycles. 


Figure 2 . 6 



AREA MJ/m 



Variation of the loop area (per cycle plastic 
work) with the plastic strain range. The loop 
areas calculated from the ELKP-parameters : U 1 = 
2* eP*av, are close to the actual values. 


Figure 2.7 
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Specimen 

ID 

eP 

T m 

(MPa) 

n f 

N 

*k 

(MPa) 

°oc 

(MPa) 

CT a 

(MPa) 

M 

(GPa) 


D7 

0.0005 

0 

608 

500 

926.0 

978.9 

1671.7 

193.5 


D8 

0.001 

0 

284 

250 

1063.7 

715.5 

1736.3 

182 . 3 

= 

D9 

0.002 

0 

305 

300 

719.9 

-297.2 

1736.9 

167.7 


DIO 

0.001 

200 

351 

350 

746.0 

225.1 

+923 . 0 

180.6 









-642.7 






Average 

864.0 



181.0 

— 


eP - Plastic strain amplitude 
Shear mean stress 
Number of cycles to failure 

Number of cycles at near-end-of-lif e beyond which the 
loop parameters are affected by the growing crack. 
Kinematic yield strength 

Cyclic, 0.02% offset yield strength (measured from 
zero stress) 

Stress amplitude 
Plastic modulus 
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The area of the hysteresis loop corresponds with the plastic 
work, the bulk of which is converted to heat. Figure 2.7 shows 
that the area of the bilinear, ELKP-loop (U‘ = 2* eP'a^) is a 
good approximation of the actual loop area for small strain 
amplitudes. The variation of the cyclic life with strain 
amplitude can be approximated by the relation: N=A( eP/2) 111 where 
A = 7.94 and m = -0.5. 

2.4 Cyclic Stress-Strain Properties of 7075-T6 Aluminum 

Results of the hysteresis loop measurements for the hardened 
aluminum alloy are summarized in Figures 2.8 - 2.15 and in Table 
2.3. Examples of the hysteresis loops are shown in Figure 2.9. 
The horizontal ticks mark the 0.02% plastic strain offset. 
Figure 2.9 illustrates that the loop shape changes little in the 
range 20 < N < 780 stress cycles. Figure 2.10 illustrates the 
loop shapes produced by cycling with a mean stress a M = -100 MPa. 
Initially, the mean stress produces a non-fully reversed loop. 
As is the case with the 440C steel, the non-reversibility, NR, 
decays with increasing number of cycles. The loop is essentially 
fully reversed after N = 330 strain cycles. This, and the 
observations that the ELKP parameters are insensitive to mean 
stress (see Table 2.3) while the conventional parameters show a 
dependence, are consistent with kinematic behavior. Figures 2.11 
and 2.12 show that none of the loop parameters, conventional and 
ELKP, are sensitive to strain amplitude in this case. 

The kinematic yield strength, plastic modulus and stress 
amplitude change very little with increasing number of stress 
cycles as shown in Figures 2.13 - 2.15. In the absence of 
significant cyclic hardening, the near-end-of-life values listed 
in Table 2.3 appear to be useful estimates of the ELKP-parameters 
appropriate for large numbers of stress cycles. Figure 2.16 
shows that the ELKP-parameters provide a reasonable description 
of the loop area. 

2.5 Discussion and Conclusions 

The ELKP-parameters derived from the hysteresis loops of the 
440C steel are very similar to those previously found for har- 
dened 1070 and 52100 bearing steel (Hahn et al., 1990). For this 
reason, the previous average values, modified slightly to account 
for a strain rate effect: = 1050 MPa and M = 188 GPa which are 
reasonable approximations for all three materials, have been used 
in the finite element calculations of rolling contact described 
in Section 3. 

The hysteresis loop measurements performed on the 7075-T6 
aluminum reveal that while the resistance of this material to 
cyclic plasticity is 1/3 to 1/2 of that displayed by 440C steel, 
the cyclic plastic behavior is kinematic, and analogous to the 
cyclic behavior of 440C. The use of aluminum as a model material 



AXJAL STRAIN (X) 


Examples of the axial stress-strain hysteresis 
loops displayed by 7075-T6 aluminum for constant 
plastic strain range, kcP = 0.001, by the N = 20 
and N = 780 cycles . 


Figure 2 . 8 


NORMAL STRESS (MPa) 


AXIAL STRAIN (X) 


Figure 2.9 


Examples of the axial stress-strain hysteresis 
loops displayed by 7075-T6 aluminum for constant 
plastic strain range cycles, AeP = 0.0025 and a 
mean stress, a m = -100 MPa by the N = 30 and N = 
330 cycles. 


, MPa 



Number of Cycles 


Figure 2.10 Variation of the kinematic yield strength of 7075- 
T6 with number of cycles for different values of 
the constant plastic strain range. The £>eP = 
0.00025 strain range test was carried out with a 
mean stress, ct™ = -100 MPa. 





100 1000 

Number of Cycles 


Figure 2.11 Variation of the plastic modulus of 7075-T6 with 

number of cycles for different values of the 
constant plastic strain range. The k eP = 
0.00025 strain range test was carried out with a 
mean stress, a m = -100 MPa. 



, MPa 


a Ac =0.025% 
□ Ac P =0.07% 
o Ac P =0.1% 


Number of Cycles 


Figure 2.12 


Variation of the stress amplitude of 7075-T6 with 
number of cycles for different values of the 
constant plastic strain range. The AeP = 
0.00025 strain range test was carried out with a 
mean stress. o« = -100 MPa. 




and M 



Figure 2.13 Variation of the kinematic yield strength, stress 

amplitude and plastic modulus of 7075-T6 with the 
plastic strain range. The value of for zero 

strain range was obtained by a linear 
extrapolation. Note that the half equivalent 

strain range corresponds with the conventional 
strain range. 
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24 



Figure 2.14 Variation of the conventional cyclic strength 
properties of 7075-T6 (measured from zero stress) 
with the plastic strain range. Values quoted for 
a strain range of &eP = 0.00025 were obtained 

with a mean stress, = -100 MPa, which produces 
different conventional strength values during the 
forward and reverse part of the cycle. Note that 
the half equivalent strain range corresponds with 
the conventional strain range. 






Table 2.3 Near-end-of-1 ife values of the hysteresis 
loop parameters for 7075-T6 Aluminum 


26 


« 






ft a 

« — 1 

VD 

CN 

CD 


£ h 

03 

03 

03 

03 


03 






co 






- s 






cd 






E 


O 

CO 

CP 



(N 

<n 


VD 


ft ft 

• 

• 

• 

• 


£ 

o 

o 

o 

o 







03 






ft 






a 

cd 

o 

VD 




- E 

cd 

CN 

in 



ft \ 

• 

• 

# 

• 

• 

ft 

O 

O 

o 

o 


£ 





l n 












ii 


CP 


CP 

o 


03 

• 

• 

• 

• 

„ — v 

g 0, 

IT) 

in 

CD 

CN 

0 

o 

LO 

in 

m 

in 

CP 

" — ■ 





03 






P 






<D 

03 

VD 

o 

VD 


> 

03 ft 

in 

cd 

VD 

r- 

03 

b £ 

"3* 











£ 

O— 






-03 

O 

cH VD 

in 

rH 


oft 

in 

in rH 

CP 

VD 


b £ 

cd 

^ CNJ 

CN 

CN 




+ + 










03 






(DO. 

CM 

cm n* 

CO 

in 


b £ 

0> 

on cd 

CN 

o 



rH 

<N + 

cn 

CN 




+ 




N 






03 






*ft 

cd 

CM 

rH 



b £ 

o 

CO 

CO 

VD 


v “-"' 


n 

CD 

CD 



O 

O 

in 

O 


£ 


cd 


CO 



in 

CD 

CD 

r- 

03 






ft 






£ 

N 






03 


O 



CN 

ECU 

o 

O 

o 

o 

CO 

b £ 


t — 1 



CD 

- — 


1 









II 








in 

in 

o 

o 

(D 

a — 

On] 

CN 

I s - 

o 

tr> 

V o\° 

o 

o 

o 

rH 

03 


* 

• 

• 

• 

P 


o 

o 

o 

o 

0) 






> 

P 





03 

CO O 






CD £ 

rH 

<N 

CD 








b 


CD 

w 

p 

w 

P 

CD 

rH 

P 

'rH 

■P 

03 

co 

4H 

0 

Q) 

p 

P 

CD 

0 

N 

4h 


a) 

E 


0 


— P 

-p 

CO 4h 

0) 

CO 

w 

CD ft 


P CD 

03 

-P P 

-P 

CO b 

03 

cn 

T3 

0 03 


P 0) 

-P 

a> e 

10 

Cs3 '— 

03 


rH 

E X 


0 -P 

0) 

P tP 


4H c 

•p 

CD 


T3 P 

0 

0) -P 

-p 

P W 


b 


W T3 

c 

03 rH 

■H X! a) 0 } 

ft 

-P E *H 

C 

CP'— >i 

O 

c 

ft CD -P -P 

co 

p *h a) 

CD -P E (0 

CO P 

(/) ‘H P 


a) m p i— i u-i d> 


T5 

CD 

0 T3 

O T3 

CO 

b 

P 

U rH u 


^3 

b 


-P 

CD -H 

o\o -P 

i — i 

•H 

03 

P *H -p 

CN 

•H 

P 

rH 


CD >i C/3 

O 

rH 

T3 

ft C XI 03 

• 

ft O 

E 

03 

E O rH 

O 

E 

E 

03 

CD 

3 *H <D 


03 



E 

C 



0 

C 


03 O 

0 

Cfl 

-H 

•rH 

P 

CD £ *H 

•H 

ca 

-P 

03 

03 rH Q) rH 

r-H 

0) 

LO 

P 

0) 

0 C U 

o 

p 

03 

•P X 

>i-H >, >,4J 

rH 

CQ 

CO 

U ^ o 

o 

CO 

ft 


i i i i i i i i 


ft o 

p E X Q) O 
b £ o o o b £ 


m 

p 

CD 

CD 

e 

ra 

P 

fti 

a 


- < 

D ft 


Loop area calculated from ELKP- 
Actual area of recorded loop 



27 


offers several new opportunities which are explored in Section 3. 
For one thing, the cyclic plasticity generated in hardened 
aluminum can be revealed metallographically , and this offers a 
way of testing elastoplastic finite element analyses of rolling 
contact. Secondly, the subsurface contact failure proceeds in 
as little as N "10 6 contacts in hardened aluminum even at 
relatively modest relative contact pressures. These may provide 
opportunities for studying the effects of residual stresses and 
material variables on the process and for testing fracture 
mechanics analyses of spall growth. 
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3. FINITE ELEMENT ANALYSES 
3 . 1 Background 

This section describes finite element calculations of the 
plasticity and residual stresses produced by repeated contact in 
the 3-ball-rod rolling contact fatigue testing machine. The 
machine is used in this study (see Sections 5.2) and by other 
investigators to evaluate material resistance to rolling contact 
failure. Detailed analyses of the test results and their in- 
terpretation must take into account the two distinct failure 
origins described in Section 1.1, as well as the cyclic 
plasticity present at those places. 

Service contact failures can be of either the subsurface of 
near-surface type. However, the work described in Section 5 
shows that the failure of 440C steel obtained with the 3-ball-rod 
tester with standard roughened balls and ground rods are of the 
near-surface type. Consequently, the calculations in this 
section are not directly applicable to the test results. The 
calculations would be relevant for tests performed with the 3- 
ball-rod machine with lapped balls and rods that produce sub- 
surface failures. 

Three dimensional, elastoplastic calculations of repeated 
contact are complex and, until recently, relatively little 
progress has been made. A brief summary of work in this field 
follows . 

A review by Johnson (1986) examines the possible mechanisms 
of failure and predicts the nature of the residual stress state. 
Ponter et al. (1985) apply the kinematical shakedown theorem to 
investigate the mode of deformation for rolling and sliding point 
contacts. The authors calculate optimal upper bounds for both the 
elastic and plastic shakedown limits for varying coefficients of 
friction and shapes of the loaded ellipse. Bower et al. (1986) 
used the above mentioned theorem to closely look at the 
conditions under which cumulative deformation occurs in the 
corner of a railhead, assuming elastic-perfectly plastic 
behavior. The study was further extended to a work-hardening 
quarter space. Hills and Sackfield (1984) studied the yield and 
shakedown states in the contact of generally curved bodies, with 
and without friction. Hills and Sackfield (1983a, 1983b, 1986) 
have done additional work treating the point contact problem 
mathematically . 

Kalker (1979) has developed a computer code for treating 
elastic 3-dimensional rolling contact with dry friction. Kannel 
and Tevaarwek (1984) presented a computer model for evaluating 
the subsurface stresses incurred during rolling-sliding contacts. 
Hardy et al. (1971) developed a finite element model of a rigid 
sphere indenting (not rolling on) an elastic-perfectly plastic 
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half -space. Equations were obtained by Hamilton and Goodman 
(1966) for the complete stress field due to a circular contact 
region carrying a 'hemispherical' Hertzian normal pressure and a 
proportionally distributed shearing traction. Chiu and Hartnett 
(1983) have presented a numerical method of solution for three- 
dimensional Hertzian contact problems involving layered solids. 
Hills and Ashelby (1982) have analyzed the residual stresses and 
their influence on the contact load bearing capacity for 3- 
dimensional rolling. Rydholm and Fredriksson (1981) devised a 
finite element model for analyzing shakedown problems in 3- 
dimensional rolling contacts for elastic-perfectly plastic and 
kinematic hardening material responses. Martin and Hay (1972) 
developed a 3-dimensional finite element model to analyze the 
yielding of a rail material, the subsequent development of 
residual stresses, and plastic flow due to a moving load. Line 
contact of two cylinders or of a cylinder and a half-space, with 
plane strain deformation, has also been studied in detail by 
Bhargava et al. (1985a, 1985b, 1986), Merwin and Johnson (1963) 
and Johnson and Jefferies (1963) . 

Ghonem and Kamath (1984) and Bhargava et al. (1986) have 
demonstrated that the cyclic stress strain properties of rail 
steel approach elastic-linear-kinematic-hardening plastic 
(ELKP) behavior after 10 3 < N < 10 4 cycles. Kinematic behavior 
can be expected to produce significantly different continuing 
plastic deformation and residual stresses. This was confirmed by 
Bhargava et al. (1988a) in a 2-dimensional elasto-plastic finite 
element study, where it was found that the residual stresses were 
50%, and the cyclic strains were an order of magnitude smaller, 
for ELKP than for EPP calculations at comparable relative peak 
pressures. McDowell and Moyar (1986) used a Mroz type kinematic 
hardening rule (which has been proven to be reliable for non- 
proportional loading) to perform an approximate numerical 
calculation for 2 -dimensional rolling-sliding contact, 
incorporating rail steel properties. The experimental work by 
Bower and Cheesewright (1988) and calculations by Bhargava et al. 
(1988b) reinforce the view that the rail-wheel contact must be 
treated as a 3-dimensional problem. 

Repeated contact loading for highly stressed mechanical 
components is accompanied by small amounts of continuing 
plasticity. The plastic deformation produces net shape changes, 
residual stresses, and may lead to the nucleation and growth of 
cracks. The effect of contact loads have been experimentally 
studied by Groom (1983) and Bower and Cheesewright (1988). 
Though there exist a number of treatments of the 3-dimensional 
problem of rolling contact, the information obtained is limited. 
Most of the analyses employ elastic (Johns and Davies (1976) , 
Hellier et al. (1986)), or elastic-perfectly plastic (Martin and 
Hay (1972), Johns and Davies (1976), Keer et al . (1986), Orkisz 
et al. (1986)) material behavior. Some of them evaluate the 
shakedown limits and provide peak values of certain normal 
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residual stresses. However, there is very little information on 
residual stresses, strains and plastic strain range distributions 
in the half space, especially in the vicinity of the contact. 
Only two of the studies deal with the stress-strain history 
(Ponter et al. (1985), Bower et al . (1986)). 

Previous analyses have reguired arbitrary simplifying 
assumptions; for example, in theoretical studies certain stresses 
were neglected, and in finite element studies the boundaries were 
assumed to be rigid. For three-dimensional contact, all six 
components of residual stress are possible, and they are 
functions of the position of the material point with respect to 
the load. In addition, out-of-plane plasticity is expected in 3- 
dimensional contact. Overall, the cyclic stress-strain variation 
is much more complex than in a 2-diroensional configuration. 

A physically more accurate treatment of the 3-dimensional 
frictionless rolling contact problem has been presented by 
Kulkarni et al. (1988, 1990a, 1990b, 1990c), by means of a finite 
element calculation. This model addresses the elastic 
displacement of the boundaries and incorporates the 
experimentally observed elastic— kinematic— hardening— plastic 
constitutive behavior. This model is used for the present 
calculations and is described below. Du et al . (1990) have 
extended Kulkarni ' s model to account for the presence of surface 
shear tractions. Initial calculations are introduced and compared 
to those resulting for the pure rolling case. These methods are 
adapted here to the 3-ball-rod tester contact geometry. 

The mechanisms of elasto-plastic rolling-plus-sliding 
contact with friction, in the absence of heating, has been 
examined by Johnson and Jefferies (1963) and Ham et al. (1988). 
In reality, the heating that accompanies friction will introduce 
thermal stresses and locally alter the elastic and plastic 
properties of the material. Existing literature predominantly 
deals with a) surface limited analyses: Barber (1971a, 1972, 
1973a, 1976, 1980a, 1980b, 1982), Comninou and Dundurs (1979), 
Comninou et al. (1981), Hills and Barber (1986), Korovchinski 
(1965), Mikic (1974), Panek and Dundurs (1979), and b) the 
contact of rough surfaces, including flash temperatures (the 
localized transient temperatures due to asperity contact) : 
Archard (1959), Blok (1937, 1963), Holm (1948), Kuhlman-Wilsdorf 
(1985), Nagaraj et al. (1979), Winer and Cheng (1980) and Zumgahr 
(1987). 

Many studies have examined the competing processes of: a) 
frictional heating at the contact interface and the resulting 
thermal expansion with b) heat transfer and wear. The unstable 
increases in the contact stress, known as thermoelastic 
instability (TEI ) have been investigated by: Barber (1967, 1968, 
1969, 1971b, 1973b), Burton et al. (1973), Dow (1972), Dow and 
Burton (1972, 1973), Dow and Stockwell (1977) and Johnson et al. 
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(1988). The manifestation of TEI as scuffing has been studied by 
Durkee (1978), Durkee and Cheng (1979), and Sovak and Cheng 
(1982). Except for Newman (1986), who has developed an elasto- 
plastic finite element model to study the stress field in welded 
plates, all the investigations discussing the subsurface effects 
of thermomechanical contact are elastic: Kennedy and Karpe 
(1982), Ling and Mow (1965) , "Mow and Cheng (1967) and Tseng and 
Burton (1982) . 

Premature failure due to overheating of HPOTP bearings used 
in the Space Shuttle Main Engine (Maurer and Pallini (1985), 
Bhat and Dolan (1983)), confirms the need for an elasto-plastic 
thermo-mechanical finite element model to simulate frictional 
rolling contact more realistically. A preliminary step in that 
direction was taken by Kulkarni et al. (1988) with the 
introduction of a stationary thermo-mechanical analysis of a 
conductive half plane with a convective surface subjected to the 
combination of mechanical pressure and different thermal loads 
(representing the heat generated due to friction) . Although the 
mechanical effect of the surface tractions and the cyclic nature 
of the applied loads was not considered, the study established 
the existence of varying levels of residual tensions proportional 
to the attending temperature gradients. The work describes a 
transient, translating, elasto-plastic, thermo-mechanical finite 
element model of a 2-dimensional frictional rolling contact. The 
temperature variations, stress-plastic strains distributions and 
deformations are calculated Tor “specific contact conditions. 
Additional work, now in progress, includes modelling of the 
three-dimensional 3-ball-rod pure rolling contact experiments, 
and three-dimensional rolling-plus-sliding studies. 

The present work reports the results of several finite 
element studies. Three-dimensional elasto-plastic finite element 
models have been developed to study frictionless rolling and 
rolling-sliding contact on a 3-ball-rod contact fatigue tester 
configuration. These analyses have been performed for an AISI 
440C steel and a 7075-T6 aluminum alloy. In addition, a two- 
dimensional finite element mesh was developed to simulate 
rolling-plus-sliding contact with heat generation. Descriptions 
of the models and results of the analyses are presented below. 

3.2 Analytical Procedures 

Frictionless rolling of a sphere on a semi-infinite body is 
simulated by translating a semi-ellipsoid pressure distribution 
over one face of a three dimensional finite element mesh. The 
finite element model is intended to reproduce the loading 
conditions prevailing on the 3-ball-rod experimental set up 
described by Glover (1982). A finite sized mesh, Figure 3.2.1, 
is used to represent a semi-infinite body by applying the 
appropriate elastic displacements on the other faces of the 3- 
dimensional mesh. For this purpose the semi-ellipsoidal pressure 


Figure 3.2.1 Finite element mesh used for the 3-dimensional 
calculations, the size of the elliptical contact 
patch and the axis (with their nomenclature) are 
also indicated in the figure. 




distribution is discretized into numerous concentrated forces 
directed along the z-axis. 
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Rolling-plus-sliding is represented by adding to the above 
mentioned pressure a distribution of shear tractions resulting 
from considering uniform slip at the contact. The shear 
tractions are obtained by first determining the "consistent nodal 
forces" resulting from the application of the semi-elliptical 
pressure distribution. These forces are the reaction forces which 
would arise at the nodes on the top surface, were these 
constrained to move along the z-direction. These nodal reactions 
are then multiplied by a constant friction coefficient and 
applied along the y-direction. 

The displacements of the boundary nodes are then calculated 
for each of these concentrated forces along the y- and z- 
directions using the appropriate solution to the Boussinesq's 
problem, and then superimposing the solutions for many loads (the 
relevant equations used for these calculations are presented in 
Appendix 1). Kulkarni et al. (1990a) have tested the accuracy of 
the model against a closed form analytical solution for the 
purely elastic indentation of a 3-dimensional space; there is 
excellent agreement between the results. The same authors have 
shown that the presence of the 'elastic 1 boundaries does not 
affect the distribution of residual stresses. 

Because of the symmetry of the pressure distribution about 
the rolling direction, the 3-dimensional mesh has the shape of a 
quadrant of a circle extended in the second dimension. The 
present study, in the most general case, pertains to two 
mechanical components in a three dimensional rolling situation, 
and hence the x-axis will be referred to as the axial direction, 
the y-axis, or the direction of translation of the load, as the 
circumferential direction, and the z-axis as the radial 
direction. The mesh is finer at the center, and coarser towards 
the 'elastic' boundaries. The mesh is lOw^ long in the 
circumf erential direction, and extends 5w^ along the radial and 
axial directions, where w^ is the dimension of the semi-major 
contact patch. 

The mesh is made of 8-noded linear brick elements closer to 
the boundaries, and more refined 20- and 27-noded quadratic 
brick elements toward the center. There is a total of 1392 
elements and 4649 nodes. The elasto-plastic material properties 
for steel AISI 440C and aluminum alloy 7075-T6, as well as the 
different loading schemes used in the calculations, are shown in 
Table 3.2.1. Table 3.2.2 shows the different loading conditions 
and geometries studied for the present report. 

Three loading levels were analyzed for the pure rolling case 
on steel, they will be hereafter referred to as high, medium and 
low pressures. For the rolling plus sliding study on steel, only 
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Table 3.2.1 Material Parameter for 3-dimensional model 


Parameter AISI 440C A1 . 7075-T6 


Young's Modulus (GPa) 

Poisson's Ratio 

Kinematic Yield Strength (MPa) 

Hardening Modulus (GPa) 

Kinematic Shear Yield Strength (MPa) 


207.0 
0.3 

1050.0 

188.0 
606.0 


70.0 

0.3 

382.0 

56.43 


220.55 
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Table 3.2.2 Loading and Geometry for 3-dimensional model 


a) Pure rolling - Aluminum 


Load 

Peak Pressure 

P^/kv 4 

Semi Contact 

Width 





Major 

Minor 

(N) 

(MPa) 



wi (m) 

w 2 (m) 

47.8 

1250.0 


5.7 

1 . 79xl0“ 4 

1.02x10 

Geometry of ball 

Rl 

= 0.0127 

m 




R2 

= 0.0127 

m 


Geometry of sample 

Rl 

' = 0.004787 m 




R2 

1 = 00 




b) Pure rolling - Steel 


Load 

Peak Pressure 

Eo/iik 1 

Semi Contact 

Width 



Major 

Minor 

(N) 

(MPa) 


W 1 (m) 

w 2 (m) 

1083 . 5 

5412 . 0 

8.93 

4.0889X10 -4 

2 . 3073xl0” 4 

438.1 

3999 . 0 

6.59 

3.0237X10 -4 

1 . 7062xl0~ 4 

96.7 

2413 . 0 

3.98 

1.8276X10" 4 

1.0312X10" 4 

c) Rolling- 

plus-Sliding - 

Steel 



1083. 5 

5412 . 0 

8.93 

4.0889X10 -4 2 

. 3073X10 -4 

Ratio of 

tangential to 

normal 

force = 0.1 



Geometry of ball R1 =0.0127 m 

R2 = 0.0127 m 


Geometry of sample 


Rl' = 0.0048 m 

R2 ' = oo 


4 p Q / k k = Ratio of Hertzian Peak Pressure to Kinematic 
Shear Yield Strength. 
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the high pressure condition was analyzed. The ratio of the shear 
traction to the normal force was assumed to be equal to 0.1. 
Only one load level was analyzed in the case of the aluminum 
alloy. 

The following assumptions are made for the analyses: a) the 
contact area, which is calculated from the elastic material 
properties and the geometry of the bodies, does not change after 
the semi-infinite body begins to deform plastically, and b) the 
applied pressure distribution remains Hertzian throughout the 
analysis. The Hertzian pressure is applied at one end of the mesh 
and incrementally translated through a distance of 8w^ to the 
other end of the mesh. The translational increments vary 
according to the position of the load on the mesh. The first 
increment, w lf is followed by two increments of 0.5w^, sixteen 
increments of 0.25W2# two increments of 0.5w^ and finally one 
increment of w^. before the mesh is unloaded. This loading 
sequence defines a single contact cycle. 

The non-linear finite element package ABAQUS (1989) was used 
for the calculations. Table 3.2.3 shows the various computational 
facilities and the required CPU times fo r each of the studies. 

The finite element mesh used for rolling-plus-sliding with 
heat generation is similar to the one used by Kulkarni et al. 
(1990a) and is shown in Figure 3.2.2. It is a 2-dimensional mesh 
consisting of 285 elements and 919 nodes. Eight noded 
isoparametric plane strain elements are used. These elements are 
biquadratic for displacements and linear for temperature 
variation. The mesh extends for 24w in the circumferential (x) 
direction and for 12w in the radial (y) direction. The mesh is 
graded; very fine in the vicinity of the global origin and 
coarser away from it. The thermo physical properties used in the 
calculations are indicated in Table 3.2.4. 

Boundary conditions 5 are applied to make the mesh simulate a 
semi-infinite half plane and to model appropriate thermal 
behavior. Side AB is the free surface; to simulate frictional 
rolling contact, a thermo-mechanical load is repeatedly 
translated across AB. The mechanical load is the combination of a 
normal Hertzian pressure distribution, and a tangential surface 
traction component. The tangential tractions are related to the 
normal pressures through the friction coefficient. The thermal 
loads are coupled to the mechanical loads (tangential tractions) 
by the velocity term. Side AB loses heat by convection; the film 
coefficient is indicated in Table 3.2.4. 

The non-surface boundaries are displaced elastically to make 


5 The derivation of the appropriate boundary conditions is 
presented in Appendix 2. 
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Table 3.2.3 Computational Requirements 


Case 


Computer 

CPU (hrs) 

3-dimensional 

Steel 

rolling 

Cray Y-MP 6 

0.5 

3-dimensional 

Aluminum 

rolling 

Cray X-MP 6 

1.5 

3-dimensional 

Steel 

rolling+ 

sliding 

Vax 3500 workstation 7 

84.0 


6 San Diego Supercomputer Center. 

7 The rolling sliding calculations are in the development 
stage; they will be implemented on the Cray Y-MP. 
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Figure 3.2.2 Finite element mesh used for the 2-dimensional 
calculations which accounted for rolling-plus- 
sliding, with the resulting heat generation. 
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Table 3.2.4 Thermo physical properties for 2 -dimensional 

calculations plus contact loading and geometry 


Parameter 


Value 


Young's Modulus 

Poisson's Ratio 

Kinematic Yield Strength 

Hardening Modulus 

Kinematic Shear Yield Strength 

Peak Hertzian pressure 

Half contact width 

Mass density 

Specific heat 

Friction coefficient 

Thermal expansion coefficient 

Film coefficient 

Thermal conductivity 


207.0 GPa 
0.3 


1050.0 MPa 

188.0 GPa 

606.0 MPa 
3.031 GPa 

0.5 x 10 -3 m 

7850.0 kg/m 3 

550.0 J/kg °C 

0.2 

1.2X10 -5 m/m °C 

150.0 W/m 2 °C 


60.0 W/m °C 
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the mesh behave as a semi-infinite half space. The non-uniform 
mechanical loads, p and t, are discretized into several 
concentratedf orces , each acting on an infinitesimal area. 

Elastic displacements due to each of these forces are 
calculated at the non-surface boundary nodes using the equations 
derived in Appendix 1, and the principle of superposition is 
applied. For the nodes along the bottom surface, DC, both X- and 
Y-displacements are prescribed while for those on the sides AD 
and BC, only displacements in the X-direction are calculated 
using the elastic solution. The following assumptions have been 
made: (1) temperature independent mechanical and thermo-physical 
properties, (2) arbitrary heat transfer coefficient h, (3) no 
heat generation due to inelastic straining, (4) the bearing 
material is assumed to be a continuum (metallurgical 
transformations accompanying the high temperatures are not 
considered) . 

3.3 Three Dimensional Rolling Contact of Bearing Steel 

Results are presented for the loaded state in the form of 
iso-contours. Given the complicated three dimensional nature of 
the problem, views along all three orthogonal axes are presented, 
i.e. the top view (along the z-axis) , the front view (along the- 
y-axis) , and the side view (along the x— axis) . The residual 
state values are presented for a set of elements located at the 
center of the mesh as a function of the depth. The stress-plastic 
strain histories in the half space are also introduced. The 
stresses are normalized with respect to k^, the strains with 
respect to kj^/G and the depth with respect to w 1 . 

As mentioned earlier, three different load levels were used 
for the pure rolling analyses. The alternating or translating 
component of the stresses acting upon the half— space, is shown in 
the form of stress contours for the three pressure distributions 
(high, medium and low) as well as for the three possible views 
(top, front and side) , for the case when the pressure is half way 
through the second contact. The schematic drawings included with 
the figures indicate the plane where the stresses are viewed. 
All figures are shown for the second pass, however results for 
the first pass were almost identical. 

Figure 3.3.1 shows the distribution of equivalent Mises 
stresses for high, medium and low pressure 8 viewed from the top 
face of the model, i.e. along the z-axis. In this figure, as well 
as all other figures representing the loaded model (i.e. all 


8 For the top and side views the contour plots represent the 
high, medium and low pressures, from top to bottom. For the front 
view the high, medium and low pressure plots are arranged from 
left to right respectively. 
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Figure 3.3.1 Contour distribution of equivalent (Mises) 

stresses on the surface of the mesh (viewed along 
the z-axis) , for high, medium and low pressures. 
The size of the contact is indicated in the 
figure. The schematic drawing shows the viewed 
slice. The values are expressed in N/m 2 (Pascals). 
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figures except those representing residual values) , the load is 
at the center of the mesh, and the contours are directly under 
the contact pressure. The elliptical nature of the contact is 
evident from the shape of the contours (the size of the 
contactwidths is indicated on the figure) . A constant ratio of 
semi-major to semi-minor contact width was used for all the 
calculations 9 . The equivalent stress is as high as 3.11 GPa at 
the center of the contact for the high pressure case (2.33 GPa 
and 1.33 GPa for the medium and low pressure cases, 
respectively). Figure 3.3.2 shows the distribution of equivalent 
Mises stresses for the side view, i.e. along the -x-axis, for the 
three pressure levels. It is interesting to note in these 
figures that the contours which indicate an equivalent stress 
above the kinematic yield strength (a^ = 1050.0 MPa) extend up to 
the surface for the high pressure case (contours 4 and higher) , 
reach close to the surface for the medium pressure case 
(contours 5 and higher), and are completely confined to the 
subsurface region for the low pressure case (contours 6 and 
higher) . The same observations may be made with reference to the 
contours shown in Figure 3.3.3. This figure shows the front view 
of the equivalent Mises stress contours for a set of elements 
located at the center of the mesh, as indicated in the schematic 
drawing. 

The translating 3-D Hertzian pressure produces highly 
complex stress distributions in the half space. This is clearly 
evident in Figures 3.3.4 through 3.3.8, which show stress 
contours for the second contact. Figure 3.3.4 shows the contact 
stress distribution (along the z direction) on the top surface. 
The contours provide evidence of the elliptical nature of the 
contact, elongated along the x-axis (perpendicular to the rolling 
direction). Figure 3.3.5 depicts the distribution of shear 
stress, u X y, near the top surface, viewed along the z-axis. The 
peak values for the antisymmetric distribution of this stress 
component are present at the edge of the contact, and range from 
plus to minus 166.0 MPa for the high pressure case, from plus to 
minus 133.0 MPa for the medium pressure case, and from -88.8 MPa 
to 77.7 MPa for the low pressure case. The most favorable 
conditions for crack propagation driven by this stress component 
would take place on a yz plane, i.e. perpendicular to the 
surface, and close to the edge of the contact. 

Figure 3.3.6 shows the distribution of circumferential 
stress, <7yy , viewed along the -x direction. As previously 
observed by Kulkarni et al. (1990a) for the case of a contact 
elongated along the rolling direction, a small region of tensile 
stress develops at the surface on either side of the semi-minor 
contact width. For the high pressure case, the magnitude of this 
tensile stress component reaches up to 333.0 MPa. However, the 


9 


The ratio w 1 /w 2 = 1.77. 
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Figure 3.3.4 Contour distribution of direct, a zz , stresses 
viewed along the z-axis (top surface) , for high, 
medium and low pressures. The size of the contact 
is indicated in the figure. The schematic drawing 
shows the observed slice. The values are 
expressed in N/m 2 (Pascals) . 
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Figure 3,3.5 Distribution of shear stresses, o X y, v ip wed along 
the z-axis. The figure indicates the size of the 
contact as well as the values for the contours. 
The stresses are indicated in N/m 2 (Pascal) . The 
schematic drawing shows the portion of the mesh 
under observation. 
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Figure 3.3.6 Distribution of circumferential, °yy> direct 

stresses under the surface. The pressure 
distribution has been translated halfway through 
one pass. The schematic drawing indicates the 
slice of material under observation. The values 
for the contours are in units of N/m 2 (Pascals) . 
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volume of material affected by the tensile stresses is not large. 
Figure 3.3.7 shows the subsurface antisymmetric distribution of 
shear stress, <7 yz . This stress component peaks at a depth of 
approximately O.JWi. The magnitude of this stress is, however, 
almost seven times as large as the previously mentioned values 
for a xy . The peak values for this stress range from plus to 
minus 1100. 0 MPa for the high pressure case, from plus to minus 

1000.0 MPa for the medium pressure case, and from plus to minus 

555.0 MPa for the low pressure case. The most favorable 
orientation for a crack to grow with this shear stress component 
for driving force would be on planes parallel to the surface. 

Figure 3.3.8 presents the contour distribution of axial 
stress, cr xx , for a slice of material at the center of the mesh. 
There is evidence of tensile axial stress at the edge of the 
contact, however, with the exception of the high pressure case 
where the stress approaches 444.0 MPa, it is negligible. In all 
cases though, the volume of material affected by the tensile 
stress is small. Figure 3.3.9 presents the distribution of shear 
stress, cj X z , under the contact. The peak value of this shear 
stress takes place at the edge of the contact, at a depth of 
approximately 0.32W2. The magnitude of the shear stress is also 
larger than a X y, but smaller than a y Z . The peak values are 1000.0 
MPa for the high pressure case, 700.0 MPa for the medium pressure 
case, and 4 3 3.0 MPa for the low pressure case. The distribution 
is also expected to be antisymmetric but with respect to the 
plane of symmetry of the mesh, i.e. y-z. 

The plastic strain increments, dePj^ j , for each translating 
increment are used to calculate the equivalent plastic strain by 
integrating them over the history of loading, based on a von 
Mises yield function and for the appropriate kinematic hardening 
flow rule. The equivalent plastic strain, eP e q, is thus a measure 
of the degree of plasticity attained at a material point. Figure 
3.3.10 shows the equivalent plastic strain contours half way 
through the second translation of the pressure distribution, 
viewed along the z-axis (top surface). Figure 3.3.11 shows the 
equivalent plastic strain distribution from a side view. The peak 
plastic activity takes place in all cases at a depth of 
approximately 0.25W]_. For the high pressure case, the peak 
equivalent plastic strain reaches values as high as 0.004; for 
the other two cases the maxima are equal to 0.00233 and 0.00044 
for the medium and low pressures, respectively. Figure 3.3.12 
shows the distribution of equivalent plastic strains viewed along 
the y-axis, it is evident that as the pressure increases, so 
does the volume of material affected by irreversible plasticity. 

Another indication of the degree of plasticity is the 
plastic work, given by the integral of the product a^j.deP^j 
over the loading history of the material point. This 
irrecoverable work leads to damage accumulation by exhaustion of 
ductility, which in turn is responsible for crack initiation and 
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Figure 3.3.10 Distribution of iso-contours for the equivalent 
plastic strains, viewed along the z-direction. The 
dimension of the contacts are indicated in the 
figure. The schematic shows the slice of material 
under observation. 
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Figure 3.3.11 Distribution of equivalent plastic strains under 

the surface, viewed along the x-direction. The 
schematic drawing shows the slice of mesh under 
consideration. 
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propagation, thus leading ultimately to failure. Figure 3.3.13 
shows the contours of plastic energy viewed along the x-axis. 
Figure 3.3.14 presents the subsurface distribution perpendicular 
to the z-axis; this figure indicates that the maximum 
dissipation of energy by irrecoverable plastic work takes place 
at a depth of 0.3 w;l in all three cases. Moreover, the maximum of 
plastic work occurs directly under the contact (centered with the 
pressure distribution) . However, the volume of material exposed 
to cyclic plasticity increases with the peak pressure. 

In the case of two dimensional rolling, line contact (a 
rigid cylinder rolling over a flat half-space) , the flat surface 
remains flat during rolling and the mechanically admissible 
residual stress state permits only two components, a r xx and a r yy, 
which vary as a function of the distance from the surface. The 
deformation mechanism is one of ratchetting parallel to the 
surface. In the case of point contact, which is a three 
dimensional problem, all six components of residual stress can 
exist, and are present. Their magnitude is not only a function of 
the depth from the surface (z distance) but also a function of 
their distance from the center plane (x distance) . The ploughing 
of the material introduces a strain gradient in the axial 
direction; this region is surrounded by material which has not 
been plastically deformed. Upon unloading, this material tries to 
recover its original shape, thus resulting in residual stresses. 

Figures 3.3.15 through 3.3.17 show the distribution of 
direct residual stresses at integration points located at a 
distance of 0 . 125w x from the yz plane, at the center of the mesh. 
The residual stresses are normalized with respect to the 
kinematic shear yield strength, k^, and are plotted as a function 
of the depth below the surface, normalized with respect to the 
semi-major contact width, W]_". Each plot presents the results for 
the three contact pressures used in the calculations. Table 3.3.1 
compares the peak direct residual stresses obtained for the 
present calculations with previous values reported in the 

literature . 

Figures 3.3.18 through 3.3.21 show iso-contours of the half 
equivalent plastic strain range for the three different loadings 
(high: p 0 /k~9, medium: p 0 /k«6, low: p 0 /k«4) for pure rolling, and 
for rolling-plus-sliding at the high load level (Fig. 3.3.21). 
The contours are for a slice of material located at the center of 
the mesh, as shown on the schematics. Figure 3.3.22 shows a 
comparison of the variation in the half plastic strain range with 
depth, for the three different loads under pure rolling, and the 
high load under rolling plus sliding (these values are taken at 
the same locations as Figs. 3.3.18-21). 
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Figure 3.3.14 Contour distribution of energy dissipated by 
plastic work, viewed along the z-direction on the 
slice of material indicated in the schematic 
drawing. The contours are in units of N.m 
(Joules) . 
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Figure 3.3.15 Distribution of residual direct stresses at the 
center of the mesh in the axial direction 
(normalized with respect to the kinematic shear 
yield strength) as a function of the depth 
(normalized with respect to the semi-major contact 
width) . Results for the three different loadings 
are indicated in the Fig. 
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Figure 3.3.16 Distribution of residual direct stresses at the 
center of the mesh in the circumferential 
direction (normalized with respect to the 
kinematic shear yield strength) as a function of 
the depth (normalized with respect to the semi- 
major contact width) . Results for the three 
different loadings are indicated in the Fig. 
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Table 3.3.1 Direct residual stresses 


Reference 

* r xx 

a r yy 

<7 r z z 

Groom (1983) 

-150 

-180 

-30 

Bower and Cheesewright (1988) 

-150 

-220 

-60 

Bhargava (1988), P 0 / k k =5 

-70 

-63 

NA 

McDowell and Moyar(1986) 

NA 

-212 

NA 

Kulkarni 10 (1990a), p 0 /k) C =9 . 2 

-115 

-180 

-40 

Present Work: 




Po/ k k =8 - 93 

-36 

-59 

-16 


14 

6 

14 

pure rolling p 0 /kj < =6.59 

-32 

-56 

-16 


19 

7 

18 

Po/ k k =3 • 98 

-16 

-32 

-7 


7 

5 

6 

roll. + slid. Po/ k k =8,93 

-45 

-121 

-68 

110 

247 

58 

Note 1: The peak compressive 

and tensile 

residual stresses are 

presented for the three different contact 

pressures 

under pure rolling, 

and also for 

the rolling 

plus 0.2 


percent sliding case. 

Note 2: All the stresses are expressed in MPa, and for the 

present work, for the second pass of the load. 


The contact ellipsoid was elongated along the rolling 
direction for this 3D calculation. 




Figure 3.3.18 Iso-contours of the half equivalent plastic strain 

range for pure rolling with high loading (po/k«9) . 
The contours are for a slice of material located 
at the ' ‘ center of the mesh, as shown on the 
schematics. Strain values are shown magnified X 
10 3 . 




Figure 3.3.19 


Iso-contours of the half equivalent plastic strain 
range for pure rolling with medium loading 
(p 0 /k«6) . The contours are for a slice of 
material located at the center of the mesh, as 
shown on the schematics. Strain values are shown 
magnified X 10 3 . 


Figure 3.3.20 Iso-contours of the half equivalent plastic strain 

range for pure rolling with low loading (p c /k«4). 
The contours are for a slice of material located 
at the center of the mesh, as shown on the 
schematics. Strain values are shown magnified X 
10 3 . 




Figure 3.3.21 Iso-contours of the half equivalent plastic strain 
range for rolling-plus-sliding with high loading 
(p 0 /k«9) . The contours are for a slice of 
material located at the center of the mesh, as 
shown on the schematics. Strain values are shown 
magnified X 10 3 . 
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3.4 Three Dimensional Rolling Contact of Hardened Aluminum 

The Mises equivalent stress distribution for the Aluminum 
properties is shown in Figure 3.4. I 11 . This figure shows the 
equivalent stress contours on two different sections of the mesh 
in the vicinity of the contact, as indicated in the schematic 
drawings. The peak value of the equivalent stress occurs at a 
depth of approximately 0.5w below the surface. 

Figure 3.4.2 illustrates the equivalent plastic strain 
contours underneath the contact. The equivalent plastic strain, 
eP eq , is obtained by calculating the equivalent plastic strain 
increment, deP eq , for each incremental pressure translation, and 
integrating it 4 over the entire loading history. It is a 
cumulative measure of plasticity in the half space. Plastic 
strain activity extends to a depth of approximately lw below the 
surface and the peak plastic strain is located at a depth of 
approximately 0.4w, as shown in Figure 3.4.2. 

The distribution of the continuing cyclic plasticity during 
the first and second contact sequences is shown in Figure 3.4.3. 
This shows the variation with depth of equivalent plastic strain 
range, AeP, which peaks at a relative depth of z/w = 0.45. The 
magnitude of eP is non-zero up to a relative depth of z/w = 1.5 
for the first contact. However, for the second contact (or 
steady state), eP extends only to a relative depth of z/w = l.Ow. 
The magnitude of peak eP remains essentially the same between 
the first and the second contacts. 

The residual stresses developed as a result of plastic 
deformation of the rim at the end of the second contact are 
illustrated in Figures 3.4.4 through 3.4.7. Contours of axial, 
a x , and circumferential, Oy, residual stress contours are shown 
in Figures 3.4.4 and 3.4.3. Residual stresses become tensile 
very close to the surface, within a relative depth of z/w = 0.2. 
Figures 3.4.6 and 3.4.7 illustrate the variation of o x and a y 
with relative depth. Peak residual stresses are compressive and 
are located at approximately z/w = 0.9. 

The cyclic shear stress-shear strain history experienced by 
the model for the two contacts is shown in Figure 3.4.8. The 
point for which the stress-strain values are obtained is located 
at a relative depth: z=0.4w. The non-fully reversed hysteresis 
loop for the first contact is fully reversed for the second 
contact. Thus the idealized ELKP material behavior allows no 
ratchetting or unidirectional accumulation of plastic strains. 

The permanent microstructural damage in 7075-T6 alloy is 


11 The distribution of individual stress and strain 
components are presented in Appendix 3. 
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Figure 3.4.1. Contours of von Mises equivalent stress on two 

different sections of the mesh (the sections are 
schematically indicated by the side of each 
figure) . The numbers on the individual contours 
represent different equivalent values of the 
contours . 





Figure 3.4.2 Equivalent plastic strain contours representing 
the plastic strain history as the pressure 
ellipsoid completes one translation. 
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Figure 3.4.4 Contours of the out-of-plane (or axial) residual 

stress, a r x . The stress values are tensile close 
to the rolling surface (as indicated by the 
contour levels 7 through 11) . 
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Figure 3.4.5 Contours of the circumferential residual stress, 

a r y. The stress values are tensile close to the 
rolling surface (as indicated by the contour 
levels 9 through 11) . 
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Normaiirea Residual Stress . 


Figure 3.4.6 Variation of the out-of-plane (or axial) residual 
stress, a r x , with normalized depth. The values of 
a r x are obtained from integration points located 
near the middle of the mesh, where the effects of 
the boundary are not significant. 
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Figure 3.4.7 Variation of the circumferential residual stress, 

a r y, with normalized depth. The values of a r y are 
obtained from integration points located near the 
middle of the mesh, where the effects od the 
boundary are not significant. 



Shear Stress, <r n |^Pa] 



In-plane shear stress-shear strain history 
experienced by a point in the model at a relative 
depth: z = 0.4 w. This is also the location of the 
maximum continuing plastic strain, eP. 


Figure 3.4.8 
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caused by continuing plasticity in the rim. The damage is 
characterized by the early formation of slip bands underneath 
the contact, resulting in a gradual destruction of the grain 
structure followed by sub-surface crack initiation. As 
illustrated in Figures 3.4.9 and 3.4.10, the shape and the extent 
of microstructural damage correlate very well with the equivalent 
plastic strain distribution obtained from finite element 
calculations . 

Sub-surface cracks occur within the plastic zone at depths: 
0.1 < z/w < 1.2 underneath the running track (here, the 
normalizing parameter is the semi-minor axis of the contact 
ellipse, w = 0.102 mm). The distribution of sub-surface cracks 
as a function of depth is shown in Figure 3.4.11. The equivalent 
plastic strain variation drawn to the same scale is also shown in 
Figure 3.4.11. While the extent to which cracks occur below the 
surface correlates well with the extent of continuing 
plasticity, the majority of the cracks do not occur at z/w = 0.4 
where the plastic strain peaks. Instead, a large number of 
cracks are found closer to the surface at z/w = 0.20. The 
distribution of crack length below the surface is shown in Figure 
3.4.12. Again, the longest cracks are located at z/w = 0.20, 
while the average length of cracks at z/w = 0.40 is approximately 
lOx smaller. A plausible explanation of the high occurrence of 
long cracks at z/w = 0.20 from the surface is that the 
compressive residual stresses decrease in magnitude at this 
depth. For instance, the circumferential and the out-of-plane 
residual stress distributions, illustrated in Figure 3.4.12 shows 
that high levels of compressive residual stresses at depths: 0.30 

< z/w < 1.2 gradually diminish and become tensile at depths: z/w 

< 0.2. Therefore, the process of crack initiation and growth 
resulting from peak plastic activity is impeded by the presence 
of high compressive residual stresses. 

The following sequence of events occurs in 7075-T6 aluminum 
before the final spall failure. Directly underneath the contact 
at depths: z w l.Ow, slip lines gradually develop as a result of 
continuing cyclic plasticity in the initial stages, i.e., N < 10 3 
cycles. These slip lines probably represent plastic activity in 
the primary slip systems. As the plastic damage continues, i.e., 
10 3 < N < 10 6 , more slip systems are activated. Gradually the 
slip bands become less discrete giving a destroyed appearance to 
the grain structure. The only visible slip bands at this stage 
are either very close to the surface: z/w < 0.1, or at the outer 
edges of the damaged region: z/w > 1.2, where the plastic strain 
tapers off. Small cracks nucleate within the damaged region: 
0.10 < z/w < 1.2. Smaller levels of compressive residual stress 
favor both the initiation and the growth of cracks at z/w ~ 0.2. 
As one of the crack branches reaches the surface, a chunk of 
material breaks off, leaving behind a spall pit. As this 
process continues, the spall pit grows larger and deeper into the 
rim. These steps are schematically illustrated in Figure 3.4.13. 



Figure 3.4.9 


Comparison of 
f rom finite 
microstructural 
after 10- 


plastic strain contours obtained 
element calculations and the 
changes produced by rolling: (a) 

contacts, and (b) after 2.3 x 10® 
contacts. The microstructure and the adjoining 
contours are shown on the same scale. 
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Figure 3.4.12 Distribution of sub-surface crack length with 
depth below the surface. The variation of the 
axial (a x ) and circumferential (ay) residual 
stresses are indicated by the dashed lines. 
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Figure 3.4.13 


RD. 

> 



C*>) 






Schematic diagram of spall formation in 7075-T6 
aluminum. 
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The implication of these results is that compressive 
residual stresses retard the process of sub-surface crack 
nucleation and growth even in the presence of continuing 
plasticity. Hence, longer lives can be obtained by artificially 
installing residual compressive stresses in a relatively shallow 
region under the surface, i.e., z/w < 0.25, or by altering the 
properties of this shallow layer. The rolling contact 
performance is strongly determined by the cyclic stress-strain 
response of the material . In spite of the microstrucfcural 
differences between bearing steels and 7075-T6 aluminum, the 
cyclic stress-strain response of these two materials is 
qualitatively similar. Therefore, the extension of the 
relationship between residual stresses, cyclic plasticity and 
sub-surface crack nucleation from 7075-T6 aluminum to bearing 
steels is justifiable. 

Some of the major limitations of the present work are 
noteworthy. Cyclic axial tests do not reproduce the non- 
proportional straining and rotation of principal shear directions 
that result from repeated rolling contact above shakedown. When 
the principal shear directions do not rotate, the cyclic damage 
is limited to a fewer number of favorably oriented slip systems 
and the damage is more intense. In addition, the absence of 
hydrostatic stresses and the larger stressed volume shorten the 
life in cyclic axial fatigue tests. This is a likely explanation 
for the 10 3 -fold difference between the rolling contact life and 
the cyclic axial tests. Also, in 7075-T6 aluminum, the 
precipitate particles are important contributors to strength 12 . 
However, the effect of particles is ignored in the current 
study . 

Future extension of the current study must treat (a) the 
strain rate effects on cyclic stress-strain response, (b) non- 
linearity of the kinematic hardening behavior, (c) the effect of 
particles, (d) thermal effects, and (e) anisotropy and texture 
effects. More work is needed to characterize the process of slip 
band appearance at a sub-microstructural scale. 

3.5 Two Dimensional Rolling-Plus-Sliding with Heat Generation 

Figure 3.5.1 shows the normalized residual stresses as a 
function of normalized depth, y/w (after unloading followed by 
cooling to the ambient temperature) . Both the axial (a z ) and 
the circumferential (a x ) stresses are found to be tensile up to 
a depth of approximately 0.5 y/w. Figures 3.5.2 and 3.5.3 
compare the circumferential and axial residual stresses 


12 Indeed, during the axial fatigue tests, a noticeable 
difference between the elastic slope in tension and compression 
was observed. This is probably an effect induced by the 
precipitate particles. 


83 


respectively for the thermo-mechanical and mechanical loading 
cases. significant differences in magnitude and nature are 
observed. Fig. 3.5.4 presents the residual equivalent plastic 
strain contours for a section of the mesh, and Fig. 3.5.5 shows 
the residual equivalent plastic strain variation with depth. 
Peak plastic activity occurs slightly below the surface. Figure 
3.5.6 shows a section of the residual state of the mesh (after 
cooling to the ambient temperature) . The surface temperatures 
half way through the passes are shown in Fig. 3.5.7. The 
temperature distribution is asymmetric about x/w = 0.0. Fig. 
3.5.8 shows the temperature contours half way through the third 
pass. Figures 3.5.9 and 3.5.10 present the hysteresis loops for 
the second and third passes, respectively. The loop is found to 
be closed after the third pass (second thermo mechanical pass) 
indicating fully reversed plasticity. 

The analyses reveal that the incorporation of ELKP material 
properties greatly reduces the plastic strains generated compared 
with those resulting from perfect plasticity. The peak residual 
equivalent plastic strain (€p) obtained here, 2.5x10 (p o /k=5.0, 
ju=0 . 1 ) , is equal to 1/25 of the value obtained by Kulkarni et 
al. (1989): 6.0xl0~ 2 (p o /k=5.0, n=0.1), for elastic perfectly 
plastic behavior under the same operating conditions. From Fig. 
3.5.5 it can be seen that heating does not have a large effect on 
the residual strains, but the peak strain occurs at a lesser 
depth. Further, it should be noted that no plastic strains 
extend beyond a depth of 2w (see Figs. 3.5.4 and 3.5.5). 

A similar trend is seen for the residual stresses. The 
normalized stresses are found to be tensile, consistent with the 
findings of Kulkarni et al. (finite element model, 1988), and 
Muro et al. (X-ray measurements, 1973) . But the magnitudes are 
far smaller than those obtained by Kulkarni et al. who reported a 
peak normalized circumferential residual stress of 1.7, clearly 
due to the difference in material properties. The residual 
stresses are found to be highly tensile and less compressive as 
compared to the non-heating case. The ratio of peak tensile 
stress to peak compressive stress is close to 2.5 for circum- 
ferential residual stresses, and 2.0 for axial residual stresses. 
A value of 2.0 has been reported by_ Muro et al. (1973), 
determined by X-ray measurements (see Figs. 3.5.1 and 3.5.2). 
The peak circumferential residual stresses are found to be about 
1.4 times the peak axial residual stresses (Fig. 3.5.3). It 
should be noted that the residual stresses vanish at a depth of 
2w as do the residual strains. The highest temperature after 
three passes is found to be 590 C. The temperature increment per 
pass decreases with each pass but the thermal steady state is not 
reached (Fig. 3.5.7). The temperature contours half way through 
the third pass show high temperatures away from the position of 
the load in the direction opposite to the translation of the 
load. Temperature gradients are seen to vanish after a depth of 
0.5w. This steep temperature gradient produces very high 
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Figure 3.5.4 Residual equivalent plastic strain contours 
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Figure 3.5.8 Temperature contours half way through the third 









94 


compressive stresses. This is followed by a nonuniform thermal 
contraction during cooling which appears to be the reason for the 
residual tensile stresses. These transient residual tensions 
(due to the translating thermo-mechanical load) suggest a 
possible mechanism for thermo-cracking. The hysteresis loops show 
that plasticity is fully reversed after two passes as reported 
by Ham et al. (1989). 

Temperature dependence of thermo-physical properties is 
complicated and difficult to implement. Of the thermal 
properties, r , the thermal diffusivity, usually co-varies with 
the thermal conductivity, K, at intermediate and high tempera- 
tures. As a rule p decreases and c p increases with increasing 
temperatures, while k increases or decreases with temperature, 
depending on the material. The most important thermal parameter 
is h ( T ) , the coefficient for forced convection (film 
coefficient) . In forced convection h(T) varies from 10 to 10 4 
W/m 2 K, and depends on geometry, flow conditions and physical 
properties. An attempt at a boundary layer analysis to solve the 
heat transfer problem and define h(T) for a specific contact 
situation is recommended. Of the mechanical properties, E, 
Young's modulus, decreases, while v, Poisson's ratio, increases 
with increase in temperature. 

Accounting for the changes in mechanical properties with 
temperature is recommended for further study; the material will 
soften at high temperatures and this could significantly alter 
the residual stresses and strains. 

3.6 Conclusions 

Three dimensional finite element calculations of repeated 
frictionless rolling contact have been carried out for ELKP 
properties of an AISI 440C bearing steel, at relative Hertzian 
pressures of p 0 /k k ~ 4.0, 6.6 and 8.9, using the geometry of the 
3-ball-rod contact fatigue testing machine. The subsurface 
stresses, plastic strains, residual stresses and plastic work 
done are evaluated. Equivalent stresses above the kinematic yield 
strength extend to the surface only for the highest load level. 
The peak equivalent plastic strain occurs at approximately 0.25wi 
(w 1 /w 2 = 1.77) for all load cases; the values of these strains 
are 4.4xl0“ 4 , 2.33xl0 -3 , and 4.0xl0“ 3 , for the low, medium and 
high loads, respectively. 

Three dimensional finite element calculations of repeated 
frictionless rolling contact have been carried out for ELKP 
properties of the 7075-T6 aluminum alloy, at a relative Hertzian 
pressure of Po/^k ~ 5.7, using the geometry of the 3-ball-rod 
contact fatigue testing machine. The subsurface stresses, plastic 
strains, and residual stresses are evaluated. The peak plastic 
strain amplitude, AeP max /2 = 1.6xl0 -3 , occurs at a relative depth 
of z/w = 0.4 below the surface. Peak residual stresses obtained 
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are a r x (max) = 26.4 MPa, and a r y(max) = 39.7 MPa, occuring at a 
relative depth of z/w = 0.9. 


Two dimensional finite element calculations of repeated 
rolling-plus-sliding contact have been carried out for ELKP 
properties of an AISI 440C bearing steel, at relative Hertzian 
pressures of Po/k^ = 5.0. Plastic strains, residual stresses and 
temperature distributions are presented. Plastic strains are 
shown to be considerably lower for ELKP properties than for an 
elastic-perfectly-plastic material (Kulkarni et al., 1989). 
Steep temperature gradients close to the surface produce high 
compressive stresses; the ensuing cooling causes residual tensile 
stresses, a possible mechanism for thermo-cracking. Thermo- 
mechanical behavior is highly dependent on the material 
properties, which also may be temperature dependent. A more 
precise definition of these properties, especially the thermal 
film coefficient, is recommended. 
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4. EVALUATION OF THE FRACTURE MECHANICS DRIVING FORCE FOR SPALL 
GROWTH 


4 . 1 Background 

This section reviews and extends calculations of the frac- 
ture mechanics driving force for cyclic crack growth under 
rolling contact (see Keer and Bryant, 1983, O'Reagan et al. f 
1985, Bastias et al., 1989, and Hanson and Keer, 1991, for 
reviews of the problem) . Observations of the spall process in 
440C samples, described in Section 5.4, as well as observations 
of spall damaged HPOTP bearings (Bhat and Dolan, 1982, 1983), 
indicate that the spalls nucleate close to the surface. The 
number of cycles needed by surface— breaking, micron-size crack 
nuclei to grow to a spall is one component of the contact life. 

The analysis of the growth component is a difficult task, 
particularly for the conditions in the HPOTP bearing. First, 
there are several features that can add to the K - crack 
growth driving force produced the Hertzian stresses. These 
include stationary pressure spikes arising from surface irregula- 
rities, the pressure of lubricant in the crack cavity, centri- 
fugal and tensile residual stresses superimposed on the Hertzian 
stresses and thermal stresses. The contributions from these 
sources are not necessarily in phase. Secondly, there are Mode 
I, II and III contributions to the driving force whose combined 
effect is difficult to predict. Thirdly, the driving force is 
modified by the crack face friction and the non-planar crack 
profile, features that are difficult to quantify. Fourth, only a 
few studies have examined the actual 3-dimensional crack con- 
figuration (Kaneta and Murakami, 1991, Hanson and Keer, 1991). 
Finally, analyses of the driving force are difficult to test 
because the centrifugal stresses, lubrication and thermal condi- 
tions in the HPOTP bearing are not reproduced in the 3-ball-rod 
testing machine. 

The sources of driving force are examined more fully in 
Section 4.2. A new analysis by Bower (1989) for a 2 dimensional 
surface breaking crack with entrapped lubricant is applied in 
Section 4.3 to define the fracture mechanics driving force for 
conditions of interest. Results of driving force calculations 
for the different sources are compared in Section 4.4. The 
corresponding threshold crack sizes for growth are examined in 
Section 4.5. Actual measurements of the growth life are 
presented in Section 5.4. Fracture mechanics predictions of 
crack growth are compared with the measurements in Section 6.3. 

4.2 Contributions to the Spall Growth Driving Force. 

4.2.1 Hertzian Stresses . Under pure rolling contact, the 
Hertzian stresses — those associated with perfectly smooth, 
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idealized contacting surfaces — are largely compressive and, by 
themselves, fail to generate a Mode I crack driving force 13 . Low 
levels of tensile stress are generated close to the surface for 
the case of rolling-plus-sliding with friction but these are not 
likely to produce significant Mode I K-values for small, as- 
nucleated cracks. Two-dimensional treatments of sub-surface 
cracks (O'Reagan et al., 1985) show that the cyclic shear 
stresses generated under rolling contact do produce Mode II and 
Mode III K-values. The peak values occur well below the surface 
at relative depths y » w (where 2w is the contact width) provided 
crack face friction is ignored. However, even the peak values 
are modest because the friction and interlocking between the 
rough crack faces impede the crack face sliding. For lubricated 
contact, the Mode II and Mode III K-values diminish as the 
surface is approached. This is also true for surface-breaking 
cracks. Keer and Bryant (1983) find that the K-values of 2- 
dimensional surface-breaking cracks are very small for crack 
lengths just beyond the nucleation stage, i.e. a » 10 jum. In 
general, the driving forces produced by 3-dimensional cracks are 
even smaller (Kaneta and Murakami, 1990). The conclusion to be 
drawn from this is that the mechanical effects of pure rolling 
contact of idealized smooth surfaces do not promote the growth of 
micron-size surface or near-surface cracks. Other contributions 
to the K-driving force must be considered to account for the 
formation of spalls. 

4.2.2 Stationary Pressure Spikes from Surface Irregularities . 
The analyses of Goglia et al., 1984, de Mul et al., 1987, and 
Elsharkawy and Hamrock, 1991 reveal that small, sub-microns-size 
irregulaities in the running track surface, i.e., asperities, 
dents and grooves, produce large, local perturbations of the 
Hertzian pressure. These are in the form of narrow, stationary 
pressure spikes that rise and fall as the contact passes over the 
asperity. The spikes can be expected to assist the nucleation of 
near-surface cracks and then intensify the Kjj- and K III" 
values generated, thereby facilitating their growth until they 
attain a size comparable to the dimensions of the surface 
irregularity. Evidence of this is contained in early work 
showing direct relations between surface roughness and contact 
life (Soda and Yamamoto, 1982) , the deleterious effects of debris 
dents (Lorosch, 1982) and the measurements reported in Sections 
5.4 and 5.5. No work has yet been done to analyze the contribu- 
tions of stationary spikes to the crack growth driving force. 

4.2.3 Entrapped Lubricant . Experiments by Way (1935), 
Dawson (1961), and Cla yton and Hill (1986) show that spalls do 


13 The Modes I, II and III correspond with the opening mode, 
in-plane mode and out-of-plane mode of crack loading, respective- 
ly. 
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not form in the absence of a lubricating fluid. The cracks 
always propagate in the direction of the load or, stated in 
another way, cracks only form with negative inclinations 
(Nakajima et al., 1983, Fujita, 1984). Also, when one rolling 
element drives another, the driven element fails first (Soda and 
Yamamoto (1982), Nakajima (1983)). Hanson and Keer (1991) have 
rationalized these observations, noting that under the above 
conditions, the mouth of a surface-breaking crack is open and can 
accept lubricant before closing under the action of the contact 
pressure. The effects of lubricant were treated by Keer and 
Bryant (1983), Kaneta et al. (1985,1987) who assumed that the 
entrapped lubricant transmits the pressure from the mouth of the 
crack to the crack tip. This mechanism produces both Mode I and 
Mode II crack driving forces. Recently, Bower (1988) has 
employed the distribute dislocation technique of Keer and Bryant 
(1983) to calculate the stress intensity factors for 2 dimen- 
sional surface breaking cracks. Bower proposes that the crack 
mouth closes as the contact moves over the crack, pushing en- 
trapped liquid into the crack and forcing the crack faces apart, 
producing a significant Mode I crack driving force. This mecha- 
nism is sensitive to the direction of rolling. The implications 
of the Bower analysis for short, surface breaking cracks are exa- 
mined in Section 4.4. A more recent treatment of 3-dimensional 
cracks with entrapped lubricant has been reported by Kaneta and 
Murakami, 1991. 

4.2.4 Thermal stresses . The localized frictional heating 
attending rolling-plus-sliding can enhance crack growth 
directly by way of the cyclic thermal stresses and indirectly 
through thermal stress-induced-plasticity and the resulting 
tensile residual stresses. In addition, the attending reductions 
in clearances can increase the contact pressure and the heat 
generation. The latter is part of the scenario proposed by Bhat 
and Dolan (1982,1983) for failures of the HPOTP bearing. Goshima 
and Keer (1990) have analyzed the contribution of thermoelastic 
contact on the crack driving force. The influence of thermal- 
mechanical contact on plasticity-induced residual stresses is 
examined in Section 3.5. 

4.2.5 Centrifugal and Residual Stresses . Centrifugal and 
residual stresses are ordinarily not cyclic in nature. However, 
when the compressive Hertzian contact stresses are superimposed 
on centrifugal (tensile) stresses or circumferential tensile 
residual stresses, Mode I K-driving forces are obtained 
(Mendelson and Ghosn, 1986). In addition, the centrifugal and 
tensile residual stresses reduce the stress normal to the crack 
faces, and, by way of the crack face friction, increase the Mode 
II- and Mode III- components of the driving force (Chen et al., 
1988) . Consistent with this, evidence of improvements in the 
contact life produced by compressive residual stresses has been 
reported by a number of workers (Fujita and Yoshida, 1978, 
Clark, 1985, Xiao et al., 1990). Elasto-plastic calculations of 



99 


rolling contact indicate that significant, circumferential, ten- 
sile residual stresses are generated when the rolling element is 
overloaded or subjected to sliding and overheated, at contact 
pressures well above shakedown (Bhargava et al., 1990, Kulkarni, 
1991) . Mendelson and Ghosn (1986) have used a boundary integral 
method to calculate the stress intensity factors for a surface 
breaking, vertical crack. Their analysis accounts for the 
centrifugal stress but assumes the crack faces to be frictionless 
and does not take the effects of entrapped lubricant into 
account. 

4.2.6 Cyclic Plasticity . Cyclic plasticity increases the 
crack tip displacements which are responsible for the cyclic 
crack growth. This effect appears to be small for hardened 
bearing steel because the plastic hardening rate is close to the 
elastic stiffness (Bastias, 1990) . However, the contribution of 
plasticity could be important at elevated temperatures and in 
regions of the bearing softened by overheating. 

4.3 Evaluation of The Mode I and Mode II Driving Force for 
Surface-Breaking Cracks with Entrapped Liquid. 

4.3.1 Analytical Procedure . As noted in Section 4.2.3., 
Bower (1989) has used the distributed dislocation technique to 
evaluate the Mode I and Mode II crack driving force for 2 dimen- 
sional, surface-breaking cracks with entrapped liquid. This 
section offers additional calculations of the driving forces for 
cracks with the usual inclinations. The calculations employed 
the techniques devised by Bower and were facilitated by software 
made available by him. The theoretical model, which is 
illustrated in Figure 4.1, consists of a 2-dimensional surface 
breaking crack of length, a, in a half space. The definition of 
the crack inclination © is that employed by Bhargava et al . 
(1990) . The Hertzian pressure distribution, and the associated 
frictional shear tractions, are applied on the surface of the 
half space, and incrementally translated from left to right. 
They are characterized by the peak Hertzian pressure, p 0 , and the 
maximum tangential traction, q 0 . Two different tractive forces 
were examined: (a) the negative or "driving" traction 
experienced by the slower rotating, "d riven " el ement or follower 
which is directed opposite to the rolling direction and (b) the 
positive or "braking" traction experienced by the faster rotating 
driver which is directed in the rolling direction. The 
calculations were performed for: (i) dry conditions, ( i i ) 
lubricated conditions with Hertzian pressure over the crack mouth 
acting in the crack cavity, ( i i i ) qo/Po = -0.01,-0.05,- 
0.10,0.01,0.05 and 0.10 (iv) crack face friction values of n c = 0 
and 0.2, (v) crack inclinations of © = -20° and -30° and (vi) 
relative crack depths 0.2 < x/h < 1.0. While the main objective 
of this work was to extend the driving force calculations to the 
threshold crack size, the Bowers algorithms proved to be unstable 
for crack depths shallower than a/w = 0.5. It should be noted 



Figure 4 . 1 



Geometry and nomenclature used for the surface 
breaking crack problem which was solved using the 
Bowers (1989) analytical solution. The load was 
translated from left to right. The traction is 
considered positive when its direction coincides 
with the rolling direction. 



101 


that negative tractions are expected to increase the crack mouth 
opening and permit more fluid to enter the crack cavity before it 
is closed by the contact pressure than is the case for positive 
tractions. However, the Bower model does not account for differ- 
ences in the crack mouth opening arising from the sign and 
magnitude of the tractions. 


4.3.2. Results . Results for dry contact are summarized in 
Figures 4.2 - 4.12. In these cases a Mode I component is not 
produced because there is no liquid to be trapped in the crack 
cavity. The calculations show that the potency of the surface- 
breaking crack, given by the normalized driving force range, 
6Kii/(p 0 -/ a ) , decreases as the crack becomes shorter and shallow- 
er. For an a=10 /im-long, n c =0.2, 9 C = -20 ° -inclined crack in a 
driven element, subjected to a Hertzian pressure, p 0 =2 . 5 GPa, a 
traction ratio, q o /p o =-0.05: AKj=0, and AKjj/ ( p Q ./a) «0 . 1 (for the 
extrapolation shown in Figure 4.11 with dashed lines. This 
means that AKjj « 0.7, which is likely to be well below the 
threshold . 

The results for lubricated contact and cracks with entrapped 
liquid are summarized in Figures 4.13 - 4.31. These show that 
the potency of the cracks increases as the crack becomes shorter 
and shallower for the a/w-range examined. The normalized driving 
force must reach a maximum with a/w and then approach zero for 
vanishingly small a/w, but the position of the maximum could not 
be defined (see Figure 4.29). While the AKj j-driving force for 
the same a = 10 )iim-long crack and contact pressure cited in the 
previous paragraph is much larger with lubricant pressure, the 
uncertainties in the interpolation (see Figure 4.29) prevent a 
numerical evaluation. The Mode I values for the lubricated case 
are about 3x larger than the Mode II values and could well be 
above the threshold. The Mode II driving force values are about 
25% greater with positive tractions than with negative tractions 
in the lubricated case, a difference that is not observed in the 
absence of lubrication. The Mode I - AK-values, presented in 
Figure 4.21, are independent of the sign of the traction and the 
values of the crack face friction. 

4.4 Comparisons of the Contributions to the Driving Force for 

Surface Breaking Cracks. 

To facilitate a comparison of the crack growth driving 
force arising from centrifugal stresses, entrapped liquid and 
thermal effects, the results of the different 2- and 3-dimensio- 
nal analyses were reduced to the relations between the normalized 
driving force and the normalized crack length and compiled in 
Table 4 . 1 and in Figure 4.32. Most of the analyses have been 
performed for crack lengths, a, comparable to the semi-contact 
width, w, i.e. a/w » 1.0, rather than the threshold crack length, 
e.g. 0.001 < a/w < 0.01. In some cases, such as small subsurface 
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Table 4.1 SUMMARY OF FRACTURE MECHANICS K DRIVING FORCE 

CALCULATIONS FOR SURFACE BREAKING CRACKS SUBJECTED TO 
ROLLING CONTACT 


Contact a _2o— M c a © Mode &K Ref. 

Conditions w p 0 b Po-/ w 


2D, Dry 
2D, Dry 
2D, Dry 
2D, Lubric. 
2D, Lubric. 
2D, Lubric . 
2D, Lubric . 
2D, Lubric . 
2D, Lubric . 
2D, Lubric . 
2D, Lubric . 
2D, Dry , Cent . 
2D, Dry, Cent. 
2D, Dry, Cent. 
2D, Dry, Cent. 
2D, Dry , TC 
2D, Dry , TC 
2D, Dry , TC 
3D, Lubric. 
3D, Lubric. 
3D, Lubric . 
3D, Lubric 
3D, Lubric 
3D, Lubric 
3D, Lubric 


0.5 

-0.01 

1.0 

-0.01 

1.5 

-0.01 

0.5 

-0.01 

1.0 

-0.01 

1.5 

-0.01 

0.5 

-0.01 

1.0 

-0.01 

1.5 

-0.01 

0.5 

-0.05 

0.5 

-0.05 

1.0 

0 

2.0 

0 

1.0 

0 

2 . 0 

0 

0.05 

0.1 

0. 10 

0.1 

0.20 

0.1 

0.2 

0 

0.1 

0 

0.2 

-0.1 

1.0 

-0.1 

0.1 

-0.1 

0.5 

-0.1 

1.0 

-0.1 


0.2 

00 

0.2 

00 

0.2 

00 

0.2 

00 

0.2 

00 

0.2 

00 

- 

00 

- 

00 

- 

00 

0.25 

00 

0.25 

00 

0 

00 

0 

00 

0 

00 

0 

00 

0 

00 

0 

00 

0 

00 

0.5 

0.5 

0.5 

2.0 

0.5 

0.5 

- 

0.5 

- 

2.0 

- 

2.0 

- 

2.0 


20 # 

II 

20° 

II 

20° 

II 

2° ° 

II 

20° e 

II 

2°; 

II 

20° 

I 

20[ 

' I 

2°' 

I 

25 I 

II 

25° 

I 

90° 

II 

90 ° 

II 

9(f 

I 

9°; 

I 

30° 

II 

30 [ 

II 

30 [ 

II 

45° 

II 

45° 

II 

45] 

I 

45° 

I 

45° 

I 

45[ 

I 

45° 

I 


0.18 

[ 1 ] 

0.33 

[1] 

0.44 

[1] 

0.54 

[1] 

0.49 

[1] 

0.57 

[1] 

1.36 

tl] 

1.36 

[1] 

1.30 

[1] 

0.68 

[2] 

1.44 

[2] 

0.97 

[3] 

1.13 

[3] 

0.09 

[3] 

0.27 

[3] 

0.078 

[4] 

0.131 

[4] 

0.209 

[4] 

0.035 

[5] 

0.053 

[5] 

0.15 

[5] 

0.25 

[5] 

0.12 

[5] 

0.20 

[5] 

0.20 

[5] 


Ref . : 


Note 1: 


Note 2: 


[1] This Report 

[2] Bower, A . F . (1989) 

[3] Mendelson, A., and Ghosh, L.J. (1986) 

[4] Goshima, T. , and Keer, L.M. (1990) 

[5] Kaneta, M , and Murakami, Y. (1990) 

The values indicated under Ref. [3] were obtained for: 
Roller Diam: 0.5 in Roller Length: 0.57322 in 


No. Rollers: 28 


Shaft Speed 

: 25, 

500 

Shaft 

- 

IR: 

2.0 

in , OR: 

2.30233 

in 

Inner 

Raceway - 

IR: 

2.3 

in , OR: 

2.6 

in 

Outer 

Raceway - 

IR: 

3.1 

in , OR: 

3.35 

in 


Interference fit of inner raceway 
onto the shaft at 0 RPM : 0.00233 in 

The values for Ref. [4] are for: thermal loading 

parameter (r) = 10, and for a Peclet number (P e ) = 100. 


TC - Thermal Contact Dry - No Lubricant in Crack 

Cent - Centrifugal Forces Lubric - Pressure in Crack 

due to lubric. 


Key 
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cracks, the relation between t he no rmalized driving force and 
crack length is parabolic 14 (O'Reagan et al., 1985) and this 
facilitates interpolations. However, this is not the case for 
the surface breaking cracks examined here and any interpolations 
to small crack lengths must be done with caution even though the 
driving force for vanishingly small cracks does approach zero. 

The results in Figure 4.32 reveal that pressure from fluid 
in the crack cavity, centrifugal stresses and thermal stresses 
each add significantly to the AK-driving force, particularly for 
relatively short cracks, a/w < 0.1. However, they are not 
necessarily additive because they do not peak at the same 
relative contact location. Fluid in the crack cavity is 
especially potent since it generates Mode I values that are even 
larger than the Mode II values. The compilation suggests that 
all but the smallest, surface breaking, lubricated cracks may 
grow in Mode I. It is also evident that the driving forces 
generated by 3 dimensional cracks are «20% of the values obtained 
for 2-dimensional cracks. Since the real spall nuclei are 3-di- 
mensional (see Section 5.4), the 3-D results of Kaneta and 
Murakami (1990) appear to be most relevant for estimating growth 
on the 3-ball-rod testing machine. The following expressions 
approximate the variation of the driving force with crack length 
for this case: 


AKj = A (a/w) B p 0 7w (4.1) 

where A =2.3 and B = 0.3. 


4.5 Evaluation of Threshold Crack Sizes. 

Insight into the significance of the driving force values 
can be obtained by installing a separate set of curves describing 
the conditions at the threshold for growth. By inserting speci- 
fic values of the peak contact pressure, p 0 , and the threshold 
value of AK, curves describing the driving force requirements 
for different threshold flaw sizes can be generated. These 
curves are overlayed in Figures 4.33-4.34 on the results of the 
driving force analyses in Figure 4.32 for 2 contact pressures and 
2 estimates of the threshold, & K thr?shold* The coincidence of 
the results from a particular driving force source and the 
threshold curves define the critical crack size for growth in 
each case. Cracks smaller than the threshold would either not be 
expected to grow, or to enlarge much more slowly as part of the 
crack nucleation process. Since crack nucleation is assisted by 
hard particles and inclusions, the critical flaw size can be 
viewed as a rough measure of the tolerable inclusion size. 


14 


AK/p 0 ya = A or AK/(p 0 yw) = A 7(a/w) where A is a constant. 



K„.2D.D,C 

K..2D.D.L 


K 11 ,2D.D,TC 


Normalized, Mode I and Mode II stress i 
factor ranges reported by different inves 
for different sources of the cyclic crac 
driving force and for different relate 
lengths. Details can be found in Table 
following abbreviations are used: 2 

dimensional, 3D - 3-dimensional, D - 
lubricated, C - centrifugal stress and T - 
stress . 
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Figure 4.33 Curves describing the conditions at the threshold 

for cyclic crack growth for either a Hertzian 
contact pressure, p Q = 2.4 GPa and a AK^hresh = 2 
MPa or p 0 = 3.6 GPa and a ARphresh = 3 MPa. 
These are overlayed on data points representing 
the normalized. Mode I and Mode II stress 
intensity factor ranges reported for different 
sources of the cyclic crack growth driving force 
in Figure 4.32. The coincidence defines the 
conditions for the onset of crack growth in each 
case. The following abbreviations are used: 2D- 

2-dimensional , 3D - 3-dimensional, D - dry, L- 
lubricated, C - centrifugal stress and T - thermal 
stress . 
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Figure 4.34 Curves describing the conditions at the threshold 

for cyclic crack growth for a Hertzian contact 
pressure, p Q = 2.4 GPa and a ARjHRESH = 5 MPa. 
These are overlayed on data points representing 
the normalized, Mode I and Mode II stress 
intensity factor ranges reported for different 
sources of the cyclic crack growth driving force 
in Figure 4.32. The coincidence defines the 
conditions for the onset of crack growth in each 
case. The following abbreviations are used: 2D- 
2-dimensional , 3D - 3-dimensional, D - dry, L- 
lubricated, C - centrifugal stress and T - thermal 
stress . 
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surfaces, a Hertzian pressure of p 0 = 2.4 GPa and a Kj -threshold 
= 5 MPa, 3-dimensional, surface breaking cracks (and inclusions) 
with dimensions, a « 10 /am can be tolerated in a perfectly smooth 
well lubricated bearing. The centrifugal stress contribution 
appears to be negligible. The defect tolerance will be reduced 
somewhat by thermal stresses arising from loss of lubrication and 
sliding. The a « 10 /xm defect dimension is larger than the 
inclusions normally encountered in bearing grade 44 0C. Consis- 
tent with this, the contact lives of 440C samples which contain 
«10 jzm- surface breaking cracks exceed N = 6*10 8 contacts when 
they are tested at p Q = 2.4 GPa in the 3— ball-rod tester. The 
cracks do not grow under these conditions. At higher contact 
pressures the same a « 10 /am surface cracks do grow into spalls 
and reduce the contact life dramatically (see Section 5.4). 

4.6 Conclusions 

Efforts have been made to extend the Bower treatment of 2- 
dimensional, surface breaking cracks with fluid in the crack 
cavity to smaller relative crack lengths, but the Bower algorithm 
proved unstable for a/w < 0.5. In the range 0.5 < a/w < 1.5, the 
values of the normalized driving force, K/p 0 /a decrease with the 
relative crack length, a/w, for lubricated contact, but increase 
for dry contact 

A compilation of K-driving force values of surface 
breaking cracks has been assembled. This reveals that fluid in 
the crack cavity, centrifugal stresses, and thermal stresses add 
significantly to the crack driving force, particularly for 
relatively short cracks, a/w < 0.1. Fluid in the crack cavity is 
especially potent, producing relatively large Mode I driving 
force values. The driving forces generated by 3-dimensional 
cracks are « 20% of the values calculated for 2 -dimensional 
cracks . 

As a result of the modest driving force generated by small, 
3-dimensional cracks, surface cracks as large as a w 10 /am would 
not be expected to grow in idealized smooth, well lubricated 
bearings with a Mode I, ^THR^SH = 5 MPa subjected to a Hertzian 
Pressure, p Q = 2.4 GPa. This is consistent with experimental 
results described in Section 5.4. 
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5. ROLLING CONTACT FAILURE 
5.1 Background 

The life of an element subjected to rolling contact loading 
is limited by the premature spalling and heavy wear of the 
bearing races and balls. The spalls are mainly surface 
initiated, proceeding in both balls and races, and propagate 
inward with distinctive geometric features. 

The spalling process is driven by the surface frictional 
tractions, Soda (1981) , and by the stationary pressure spikes 
originated at surface irregularities, Hamrock (1990) . Additional 
complications may be caused by: a) the presence of subsurface 
secondary particles or inclusions, which may lead to subsurface 
initiated cracks, b) strain induced phase transformations, i.e. 
retained austenite decomposition, carbide dissolution, c) 
lubrication conditions at the surface, or d) thermally induced 
stress fluctuations (as treated in Section 3.5), which eventually 
contribute to the Modes I and II driving forces (as indicated in 
Section 4.4). 

This section describes experiments conducted on 440C bearing 
steel in order to gain a better understanding of the mechanisms 
leading to spall failure under different conditions. 


5.2 Experimental Materials and Procedures 

The rolling experiments were carried out on the 3-ball-rod 
rolling contact fatigue tester developed by the ball and roller 
bearing group at Federal-Mogul, and furnished by NTN Bower. The 
design and operation of the tester machine has been described by 
Glover (1982) . 

The test material is AISI 440C, martensitic stainless steel, 
heat treated to a hardness of 59-61 HRc. The test samples were 
subjected to the following heat treatment: 

a) Austenitized at 1930130 S F for 1 hour at temperature. 

b) Hardened by quenching in liquid nitrogen. 

c) Tempered one hour (minimum) at 325125 B F. 

d) Cooled in air to 70110 B F. 

e) Cold soaked in liquid nitrogen for 30 minutes. 

f) Tempered for 1 hour at 325125 fi F. 

The geometry of the samples is shown in Figure 5.1. The test 



0_ .000025 
0 .00005 


Figure 5 



. 1 Test specimen design according to Glover (1982), all 
the dimensions are in inches. 
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samples were ground to a surface finish of 6.5 AA (0.165 /xm) / 
with a roundness of 1.27-1.54 jum (a more detailed study of the 
surface finish is included in Figure 5.2). Each rod can be used 
to perform up to 8 individual tests. 

The balls used in the RCF tester are made of AISI 52100 
bearing grade steel. Two different sets of balls were used during 
the tests. The first set consisted of balls with smooth surface, 
i.e. 0.34 AA (8.6xl0 -3 /u m) , and the second set consisted of balls 
which are abrasive blasted to increase the surface roughness, 
i.e. 4.28 AA (0.11 [im) . A more detailed surface analysis of both 
types of balls can be seen in Figures 5.3a and b. The chemical 
compositions for both the 440C and the 52100 steels are shown in 
Table 5.1. 

5.2.1 Rolling Contact Testing; A total of 21 individual tests 
(see Table 5.2), distributed among three specimens were 
conducted. After each test the balls were replaced to insure 
uniform testing conditions. Exxon grade 2380 turbine oil was used 
as a lubricant in all the tests. The oil was applied by drip feed 
at a rate of 10 drops per minute. A spring load of 253.6 N was 
used in all the tests. The peak Hertzian pressure was calculated 
to be equal to 5.4 GPa . 

Six tests were conducted under similar conditions to 
determine the amount of scatter in the results (Test Series 1) . 
An attempt was also made to vary the content of retained 
austenite in the specimens by soaking in liquid nitrogen for 
different lengths of time; a total of 10 tests were devoted to 
this purpose (Test Series 2) . Three more tests were conducted 
with smooth balls in order to register any differences in the 
failure life, location and size of spalls, or the surface 
appearance after the test (Test Series 3) . Two additional tests 
were conducted using rough balls in order to study the contacting 
surface (Test Series 4) . The tests were interrupted after a short 
time after starting to examine and photograph the surface under 
the scanning electron microscope. After examination, the 
specimens were inserted back into the RCF tester. The tests were 
continued until spalling failure occurred. 

The fatigue test results summarized in Table 5.2 have been 
analyzed using the two-parameter Weibull statistics. The Weibull 
estimates for the two cases are summarized in Tables 5.3 and 5.4, 
respectively. Graphical representation of the test data is 
illustrated in Figure 5.4. The procedure adopted to obtain the 
Weibull parameters is similar to the graphical method described 
by Lipson and Sheth (1973) . A computer program was written to 
evaluate the L10 and L50, as well as the characteristic lives, 
and the upper and lower bounds for the 75% confidence levels, by 
using a binomial distribution method. 
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Table 5.1 Composition of 440C and 52100 


Steel 

C 

Mn 

Si 

Cr 

Ni 

Mo 

Cu 

A1 P 

S 











440C 

1.05 

.44 

.44 

17.11 

. 10 

.48 

.04 

.02 <.022 

<.001 

52100 

1.02 

.34 

.32 

1.45 

0.13 

.05 

.16 

<.09 

<.014 


Balance Fe. 
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3.b Surface roughness analysis of a 52100 RCF ball. 

The ball has been lapped to a surface finish of 
0.34 /iin, for a peak-to-valley of 6.34 /i in. 
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Table 5.2 Result of Rolling Tests 


Test 

Series 

_1 : 

Study of the contacting surfaces and 

scatter in 




the experimental 

data 

(combined with 

Test Series 




for the 

statistical analysis 

shown in 

Table 5.3) 

Test 

Specimen 

Track 

Machine 

Hours 

Ball 

Hours to Stress 

# 

ID 


# 

Head 

in LN 2 

Lot 

Failure CyclesxlO 

1 

A1 


LI 

1 

0 

R 

22.0 

11.33 

2 

A2 


LI 

1 

0 

R 

10.4 

5.36 

3 

A 2 


L3 

1 

0 

R 

6.4 

3.30 

4 

A2 


L5 

1 

0 

R 

2.4 

1.24 

5 

A2 


L7 

1 

0 

R 

2.4 

1.24 

6 

A2 


L9 

1 

0 

R 

10.4 

5.36 

Test 

Series 

_2. : 

Retained austenite study. 



7 

A1 


L3 

1 

1 

R 

27.8 

14 . 32 

8 

hi 


L5 

1 

2 

R 

7.5 

3.86 

9 

hi 


L6 

2 

3 

R 

22.1 

11.38 

10 

hi 


L7 

1 

3 

R 

15.2 

7.82 

11 

hi 


L8 

2 

3 

R 

14.2 

7.31 

12 

A2 


L2 

2 

2 

R 

3.9 

2 . 00 

13 

A2 


L4 

2 

2 

R 

6.4 

3.30 

14 

A2 


L6 

2 

2 

R 

1.0 

0.52 

15 

A2 


L8 

2 

2 

R 

16.6 

8.54 

16 

A2 


L10 

2 

2 

R 

3.4 

1.75 

Test 

Series 

_3 : 

Initial 

study of 

ball 

roughness effect. 

17 

A1 


L2 

2 

3 

S 

30.4 

15.65 

18 

A1 


L9 

1 

3 

S 

14 . 0 

7.20 

19 

A1 


L10 

2 

3 

S 

63 . 6 

32.75 

Test 

Series 

_4 : 

Study of the contacting surfaces and 

scatter in 




the experimental 

data 

(combined with 

Test Series 




for the 

statistical analysis 

shown in 

Table 5.3) 

20 

A6 


L2 

2 

0 

R 

19.6 

10.09 

21 

A6 


L3 

2 

0 

R 

5.3 

2.7 


15 


R=Rough Lot, S=Smooth Lot 
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Table 5.3 Two-Parameter Weibull Estimate for Case 1 

Study of the contacting surfaces and scatter in the experimental 
data are combined here. The results shown under Test Series 1 and 
4 in Table 5.2 


Item 

Cycles 

Median Ranks 

1 

0. 124xl0 7 

0.08300 

2 

0. 124X10 7 

0.20210 

3 

0. 273X10 7 

0.32130 

4 

0.33 OxlO 7 

0.44040 

5 

0.536X10 7 

0.55960 

6 

0. 536X10 7 

0.67870 

7 

0. 101x10 s 

0.79790 

8 

0.113X10 8 

0.91700 

Slope 

= 1.24469 

L10 = 0.94748x10 

Intercept 

= -19.37923 

L50 = 0.43041x10 


Characteristic Life = 0.57778xl0 7 
Correlation = 0.9633787 


Table 5.4 Two-Parameter Weibull Estimate for Case 2 

Study of the effect of retained austenite in the fatigue life 
using the results shown as Test Series 2 in Table 5.2. 


Item 


Cycles 


Median Ranks 


1 

2 

3 

4 

5 

6 


0. 520xl0 6 
0. 175X10 7 
0 . 2 OOxlO 7 
0. 330X10 7 
0.386X10 7 
0.854X10 7 


0.10910 

0.26550 

0.42180 

0.57820 

0.73450 

0.89090 


Slope = 1.11034 L10 = 0.49944xl0 6 

Intercept = -16.81940 L50 = 0.27247X10 7 


Characteristic Life = 0.37904X10 7 
Correlation = 0.9825556 



WEIBULL ESTIMATE FOR VALUES IN TABLES 5.3 AND 5.4 



Figure 5.4 Weibull plot for the RCF test results presented in 

Tables 5.3 and 5.4, peak pressure, p 0 =5.4 GPa. 
Case 1 (circles) represent the statistical 
analysis of the results of the Test Series 1 and 4 
from Tables 5.2. Case 2 (triangles) are the 
results of the study of the influence of retained 
austenite on the fatigue life, i.e. Test Series 3 
from Table 5.2. 
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In order to perforin the statistical analysis of the data 
included in Table 5.2, the samples were grouped in two different 
categories: a) Tests Series 1 and 4 were combined since they 
represent similar test conditions, b) Six tests from Test Series 
2, each with two hours of liquid nitrogen immersion, were 
selected to represent similar test conditions. The results show 
that the Weibull slopes for the two cases are within the range: 
0.8 < e < 1.5 as expected for bearing steels, according to Harris 
(1984) . While the L10 lives for the two cases are comparable, the 
L50 life for Case 1 is approximately 60% higher than the L50 life 
for Case 2. This difference may be due to the transformation of 
the retained austenite after the immersion in liquid nitrogen. In 
the case of Test Series 3, the use of smooth balls resulted in a 
relatively longer failure life. 

5.2.2 Artificial Defects: Following a procedure similar to the 
one adopted by Lorosch (1982), the running tracks of 440C 
specimens were artificially indented by a 120° sphero-conical 
Rockwell "Brale" penetrator. The load applied to obtain the 
hardness indents was approximately 10 kg, and the diameter of the 
indentation was approximately equal to 0.1 mm. Spalls 
consistently initiated close to the hardness impression in all 
cases. The specimen surface near the hardness impression was 
periodically studied and photographed using the scanning electron 
microscope. The general features of the spall produced by the 
hardness indent are very similar to the 'classical' spall 
indicated by Kumar (1987) in early reports. After the crack had 
initiated, the specimen was inserted back in the RCF tester and 
run to failure. 

A total of 30 experiments were run at three different peak 
Hertzian pressures, p D = 3.3, 4.0 and 5.45 GPa. For each pressure, 
two sets of five samples were run, one without and one with a 
hardness impression. All these experiments were run with rough 
balls. Summaries of the rolling tests for the three pressure 
levels are presented in Tables 5.5, 5.7 and 5.9. In addition to 
these, Tables 5.6, 5.8 and 5.10 present the Weibull parameters 
for the tests without and with hardness impressions for the three 
cases. The results were also plotted on Weibull charts-^; Figure 
5.5 shows the results for the tests run under a peak pressure of 
p Q = 3.3 GPa, and Figures 5.6 and 5. 7. a present the results on 
the same probability paper for the tests conducted at p c = 4.0 
and 5.45 GPa, respectively. 

In order to ascertain the influence of the microasperities 
on the total life of the sample, an additional set of 13 
experiments was run at a peak Hertzian pressure of p Q = 5.41 GPa 


16 The percent of failure is expressed as a function of the 
number of stress cycles. 
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— — 

1 r.. 

(p 0 

= 3300 MPa) 

. .... .. 





Test 

# 

Specimen 

ID 

Track 

ID 

Machine 

Head 

Hardness 

Indent 

Ball 

Loti 7 

Hours to 
Failure 

Stress 

Cycles 

1 

A9 

LI 

1 

No 

R 

87.7 

45.17 x 10 6 

2 

A9 

L2 

1 

No 

R 

57.3 

29.51 X 10 6 

3 

A9 

L3 

1 

No 

R 

25.4 

13.08 X 10 6 

4 

A9 

L4 

1 

No 

R 

69.1 

35.59 X 10 6 

5 

A9 

L5 

1 

No 

R 

57.8 

29.77 X 10 6 

Table 5.5(b) Rollina exoeriments 

with hardness imoression 




i£c 

, = 3300 MPa) 



Test 

Specimen 

Track 

Machine 

Hardness 

Ball 

Spall Growth Life 

# 

ID 

ID 

Head 

Indent 

Lot 

in Stress Cyclesl 8 

1 

A9 

L6 

1 

Yes 

R 

10.15 

X 10 6 

2 

A9 

L7 

1 

Yes 

R 

7.11 

X 10 6 

3 

A9 

L8 

1 

Yes 

R 

16.07 

X 10 6 

4 

A9 

L9 

1 

Yes 

R 

9.06 

x 10 6 

5 

A9 

L10 

1 

Yes 

R 

12.52 

x 10 6 


17 R = Rough lot, S = Smooth lot. 

18 Spall Growth Life = (Total failure life with hardness 
impression) - (Life at first observation of cracks) . 


Table 5.6(a) Weibull parameters for rolling test without 

hardness impression Cp 0 = 3300 MPa^ 
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Item Cycles 


1 

0.1308 

X 

00 

o 

iH 

2 

0.2951 

X 

H* 

O 

00 

3 

0.2977 

X 

10 8 

4 

0.3559 

X 

M 

O 

00 

5 

0.4517 

X 

00 

O 

H 

Slope 

= 2 

. 11496 


Intercept = -36.78784 


Median Ranks 
0.1294 
0.3147 
0.5000 
0.6853 
0.8706 

L 10 = 0.12361 x 10 
L 50 = 0.30123 X 10 


Characteristic Life = 0.35822 x 10 8 
Mean Life = 0.31719 x 10 8 
Correlation = 0.9503689 


Table 5.6(b) Weibull parameters for rolling te st w ith 


Item 


hardness impression (p Q = 3300 MPa) 
Cycles Median Panks 


l 

0.7110 

X 

10 7 

2 

0.9060 

X 

10 7 

3 

0.1015 

X 

f— * 
O 
00 

4 

0.1252 

X 

10 8 

5 

0.1607 

X 

10 8 

Slope 

Intercept 

= 3. 

= -53. 

28743 

65850 


Characteristic Life = 
Mean Life = 
Correlation = 


0.1294 

0.3147 

0.5000 

0.6853 

0.8706 

L 10 = 0.6186 X 10 7 
L 50 = 0.1097 X 10 8 

0.12266 X 10 8 
0.11038 X 10 8 
0.9842627 


00 00 
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Table 5.7(a) Rolling experiments without hardness impression 

fp 0 = 4000 MPa') 



Test 

# 

Specimen 

ID 

Track 

ID 

Machine 

Head 

Hardness 

Indent 

Ball 

Lot 

Hours to 
Failure 

Stress 

Cycles 


1 

A8 

L4 

1 

No 

R 

19.0 

9.785 x 10 6 


2 

A8 

L5 

1 

No 

R 

48.5 

24.980 X 10 6 


3 

A8 

L6 

1 

No 

R 

21.4 

11.020 X 10 6 

ei 

4 

A8 

L7 

1 

No 

R 

25.7 

13.240 X 10 6 


5 

A8 

L8 

1 

No 

R 

30.5 

15.710 X 10 6 



Table 5.7 (b) Rolling experiment 

with hardness impression 

— 




lEC 

i = 4000 MPa) 




Test 

Specimen 

Track 

Machine 

Hardness 

Ball 

Spall Growth Life 

— 

# 

ID 

ID 

Head 

Indent 

Lot 

in Stress cycles 


1 

A8 

L9 

1 

Yes 

R 

5.974 

X 10 6 


2 

A8 

L10 

1 

Yes 

R 

7.828 

X 10 6 

— 

3 

A8 

Lll 

1 

Yes 

R 

4.017 

x 10 6 


4 

A8 

L12 

1 

Yes 

R 

3.966 

X 10 6 

E? 










5 

A8 

L13 

1 

Yes 

R 

4 . 172 

x 10 6 


& 
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Table 5. 8 fa) Weibull parameters for rolling test without 

hardness impression (p Q = 4000 MPa) 

Median Ranks 


Item Cycles 

1 0.9785 x 10 7 

2 0.1102 x 10 8 

3 0.1324 X 10 8 

4 0.1571 X 10 8 

5 0.2498 X 10 8 

Slope = 2.63138 

Intercept = -36.78784 

Characteristic Life = 
Mean Life = 

Correlation = 


0.1294 

0.3147 

0.5000 

0.6853 

0.8706 

L 10 = 0.72331 X 10 7 
L 50 = 0.14799 X 10 8 

0.17011 X 10 8 
0.15142 X 10 8 
0.9281409 


Table 5.8(b) Weibull parameters for rolling test with 
hardness impression fp 0 = 4000 MPa) 


Item 


Cycles 


Median Ranks 


1 

0.3966 x 

10 7 

0.1294 


2 

0.4917 x 

10 7 

0.3147 


3 

0.4172 x 

10 7 

0.5000 


4 

0.5974 x 

10 7 

0.6853 


5 

0.7828 x 

10 7 

0.8706 


Slope 

= 2.94662 

L 10 = 0.2749 

X 10 

Intercept 

= -45.93936 

L^q = 0.5210 

x 10 

Characteristic Life 

= 0.5901 X 

10 7 


Mean Life 


= 0.5279 X 

10 7 


Correlation 

= 0.8593127 



7 

7 
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WEIBULL VALUES FOR TABLES 5.8(A) AND 5.8(B) 
PRESSURE - Po = 4000 MPa 


X 

0 

F 

C 

F D 
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L X 
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E 


33 


3Z 
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60 
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UX 75 

0 No Hard. Impress 
LX 75 
UX 75 

A Hard. Impress. 



1000000 


LIFE (STRESS CYCLES) 


Figure 5.6 Weibull plot for the RCF test results presented in 

Tables 5. 8. a (circles) and 5.8.b (triangles), peak 
pressure, p o =4.0 GPa. The two outher lines 
represent the limits of the 75% certainty band. 
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Table 5.9(a) Rolling experiment without hardness impression 

(p c = 5450 MPa 1 


Test 

# 

Specimen 

ID 

Track 

ID 

Machine 

Head 

Hardness 

Indent 

Ball 

Lot 

Hours to 
Failure 

Stress 

Cycles 

1 

A 7 

L4 

1 

No 

R 

8.1 

4.17 x 10 6 

2 

A7 

L5 

1 

No 

R 

8.5 

4.38 x 10 6 

3 

A7 

L8 

1 

No 

R 

7.2 

3.71 x 10 6 

4 

hi 

L9 

1 

No 

R 

14.2 

7.31 X 10 6 

5 

hi 

L10 

1 

No 

R 

12.2 

6.28 X 10 6 


Table 5.9(b) Rolling experiment 

with hardness impression 




Xe< 

3 = 5450 1 

MPa) 



Test 

Specimen 

Track 

Machine 

Hardness 

Ball 

Spall Growth Life 

# 

ID 

ID 

Head 

Indent 

Lot 

in Stress 

Cycles 

1 

A7 

Lll 

1 

Yes 

R 

0.695 x 

10 6 

2 

A7 

L12 

1 

Yes 

R 

1.210 x 

10 6 

3 

A8 

LI 

1 

Yes 

R 

1.416 x 

10 6 

4 

A8 

L2 

1 

Yes 

R 

2.678 X 

10 6 

5 

A8 

L3 

1 

Yes 

R 

4.532 X 

10 6 
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Table 5 . 10 (al Weibull parameters for rolling test without 

hardness impression ( p q = 5450 MPa) 


Item 

1 

2 

3 

4 

5 


Cycles 

0.3708 x 10 7 
0.4172 x 10 7 
0.4378 x 10 7 
0.6283 x 10 7 
0.7313 X 10 7 


Slope = 3.32260 

Intercept = -51.73905 


Median Ranks 
0.1294 
0.3147 
0.5000 
0.6853 
0.8706 

L 10 = 0.29419 X 10 7 
L 50 = 0.51865 X 10 7 


Characteristic Life = 0.57913 x 10 7 
Mean Life = 0.52146 x 10 7 
Correlation = 0.928257 


Table 5.10(b) Weibull parameters for rolling test with 
hardness impression (P D = 5450 MPa) 

Median Ranks 


Item 

Cycl 

es 


1 

0.6953 

X 

10 6 

2 

0.1210 

X 

10 7 

3 

0.1416 

X 

10 7 

4 

0.2678 

X 

10 7 

5 

0.4532 

X 

10 7 

Slope 

Intercept 

= 1 . 
= -20. 

38991 

43612 


Characteristic Life = 
Mean Life = 
Correlation = 


0.1294 

0.3147 

0.5000 

0.6853 

0.8706 

L 10 = 0.4812 X 10 
L^q = 0.1866 X 10 

0.24294 X 10 7 
0.22133 X 10 7 
0.9734292 


vo 
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WEIBULL VALUES FOR TABLES 5.11 AND 5.12 
PRESSURE - Po = 5410 MPa 



LHiEND 

LX 75 
UX 75 

N.H.I. + S.B. 
LX 75 
UX 75 

H.I. + S.B. 


LIFE (STRESS CYCLES) 


Figure 5. 7. a Weibull plot for the RCF test results presented in 
Tables 5. 10. a (circles) and 5.10.b (triangles), 
peak pressure, p 0 =5.4 5 GPa. The two outher lines 
represent the limits of the 75% certainty band. 
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with and without hardness impression 19 , but employing smooth 
balls. Tables 5 .11 a nd 5. 12 show the rolling test results and the 
significant Weibull distribution parameters. The results are 
presented in Weibull charts in Figure 5.7.b. 


5.2.3 Spall Replication and Sectioning Techniques: Before an 
accurate 3 -dimensional finite element model of a crack under 
rolling contact loading can be implemented, more needs to be 
known about the geometry of the spalls. Typical spall angles and 
dimensions, as well as their relationship to the applied load, 
the geometry of the system, the lubrication, etc. , have to be 
accurately defined. Spalls created during the experimental phase 
of this program were studied using two different methods: a) a 
replication technique, and b) a mechanical sectioning followed by 
scanning microscopy observation. 

The first procedure consists of the following steps: 1) the 
sample was degreased using methanol in an ultrasonic cleaner, 2) 
a medium viscosity plastic material (vinyl polysiloxane) was then 
prepared and applied on a glass microslide, 3) the sample, with 
the spall under study facing down, was pressed against the 
plastic and held in position for approximately 15 minutes, 4) the 
sample was removed from the plastic. In order to check the 
degree of fidelity with which this method reproduces the 
topological features of the spall, the replicas were gold coated 
and observed under the scanning electron microscope. Figure 5.B 
shows photographs of a spall and the resulting gold coated 
replica . 

The samples were mounted on the adjustable stage of a Sloan 
Dektak Surface Profile measuring device 20 , which allowed up to 
100 /xm displacement in the radial direction 21 . The stylus of the 
profilometer was positioned on the replica at a distance of 
approximately 0.3 mm from the spall (this position was taken as 
the origin in the y direction) . The stage was moved along the 
axial direction at a rate of 0.01 cm/min. Figure 5.9 shows traces 
obtained with the profilometer, as well as their relative 
location on the spall. A total of 18 traces were taken for each 
spall at intervals of approximately 0.05 mm each. The traces were 
then digitized and used to create a 3-dimensional drawing of the 


19 Two samples were not used for the calculations because 
they presented lives which were completely different from the 
rest of the population. 

20 This equipment has a sensitivity ranging from 25 Angstrom 
units to a maximum of 1x10 s Angstrom. 

21 Relative to the spall, x is the axial direction, y is 
the circumferential direction, and z is the radial direction. 



Table 5.11 Rolling test results without hardness 

impression and smooth balls fp o =5410 MPa) 
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Test 

# 

Specimen 

ID 

Track 

ID 

Machine 

Head 

Hardness 

Indent 

Ball 

Lot 

Total : 
Stress 

Life in 
Cycles 

1 

440-1 

10 

2 

No 

S 

30.2 

X 

10 6 

2 

440-1 

11 

2 

No 

S 

70.3 

X 

10 6 

3 

440-1 

12 

2 

No 

S 

4.49 

X 

10 622 

4 

440-1 

5 

2 

No 

S 

89.1 

X 

10 6 

5 

440-1 

6 

2 

No 

S 

34.7 

X 

10 6 

6 

440-1 

10 

2 

No 

S 

41.6 

X 

10 6 


Weibull Distribution Parameter 

L 10 = 21.0 x 10 6 

L 50 = 51.7 x 10 6 

Mean Life = 54.5 x 10 6 

Slope = 2.09 

Correlation = 0.9386 


Table 5.12 Rolling Test Results With Hardness 

Impression and smooth balls fp o =5410 MPa) 


Item 

# 

Specimen 

ID 

Track 

ID 

Machine 

Head 

Hardness 

Indent 

Ball 

Lot 

Growth 

Stress 

Life in 
Cycles 

1 

440-L1 

5 

2 

Yes 

S 

0.57 

X 

10 6 

2 

440-L2 

4 

2 

Yes 

S 

0.54 

X 

10 6 

3 

440-L3 

4 

2 

Yes 

S 

0.67 

X 

10 6 

4 

440-L3 

5 

2 

Yes 

S 

0.26 

X 

10 6 

5 

440-L3 

8 

2 

Yes 

S 

17.9 

X 

10 623 

6 

440-L3 

9 

2 

Yes 

S 

1.24 

X 

10 6 

7 

440-L3 

10 

2 

Yes 

S 

1.19 

X 

10 6 


Weibull Distribution 

Parameter 

L^q = 2.45 X 

10 5 

l 50 = 7.07 x 

10 5 

Mean Life = 7.71 x 

10 5 

Slope = 1.7789 


Correlation = 0.9700 



22 This sample was not used for the Weibull distribution 
calculations . 

23 This sample was not used for the calculations because the 
total life is significantly different than the rest of the 
population. 
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WEIBULL VALUES FOR TABLES 5.10(A) AND 5.10(B) 
PRESSURE - Po = 5450 MPa 


r I FT.FK in 



LIFE (STRESS CYCLES) 


Figure 5.7.b Weibull plot for the RCF test results presented in 
Tables 5.11 (circles) and 5.12 (triangles), peak 
pressure, p 0 =5.41 GPa. The two outher lines 
represent the limits of the 75% certainty band. 
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Figure 5.9 Traces of the spall's surface obtained with the 

Dektak Surface Prof ilometer . The scales are indicated 
in the figure, as well as the relative position of 
the traces in the spall. 
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spall; Figure 5.10 shows this drawing as well as the 
corresponding iso-contour map 24 . 

Some effort has been devoted in the past to characterizing 
the geometry of the spall failures. Shao et al . (1987) and Soda 
(1981) have reported cracks running at 20° to 30° to the surface. 
Bhat and Dolan (1982, 1983) have detected crack inclinations of 
10° to 20°; Voskamp (1988) has indicated the presence of cracks 
running parallel to the surface. Kumar et al. (1987) , among 
others, have pointed out the following general features of the 
spalls: a) the V-shaped spalls point in the rolling (forward) 
direction, b) the cracks grow laterally and circumferentially 
making a small angle with the running surface, c) the cracks grow 
deeper as they extend in the direction opposite to the rolling 
direction, and present a stepped like character, d) the cracks 
break through to the surface at the shallow end, thus producing a 
cantilever flap supported by an "unbroken" ligament at the deeper 
end, e) due to secondary cracks running perpendicular to the main 
crack front, portions of the flap behind the crack front break 
off periodically producing wear fragments. 

In order to systematically check these observations, spalls 
produced for the present project were sectioned and observed in 
a transmission electron microscope. The sectioning was done with 
a high speed silicon carbide cutting blade; water was used as a 
lubricant. Spalls were sectioned along the axial and 
circumferential directions. After sectioning, the spalls were 
mounted in conducting resins (Conductomet) and polished with 
polishing paper grits 240, 320, 400 and 600. Pictures were taken 
in a Hitachi scanning electron microscope. 

Table 5.13 summarizes the results obtained from the 
observations 25 . Figure 5.11 shows a schematic representation of a 
typical spall and the nomenclature used for the different 
dimensions and angles. Figures 5.12.a-d show spalls viewed along 
the radial direction and, after the sectioning, in the 
circumferential direction. In the radial view, the presence of a 
micropit right before the V-apex of the spall is evident (Figure 
5. 12. a). Figures 5.12.b and c show spalls formed in the vicinity 
of hardness impressions. The top views also indicate the presence 
of a surface breaking radial crack; this can also be seen in the 
circumferential cut. The linking of this radial crack with the 
horizontal subsurface crack results in the detachment of another 


24 The 3-dimensional drawing and the iso-contour maps are 
created by a Fortran program which makes use of the subroutine 
DISSPLA from Integrated Software Systems Corp. 

25 All these results are for a peak Hertzian pressure of p G = 
5.45 GPa , for which the semi-minor and -major contact widths are 
w 2 = 0.2 2 mm and W]_ = 0.40 mm respectively. 


radial \m*\] 


0 . 2 * 


Sample A 7 
Spall L5 



Figure 5.10 3-Dimensional digitized view of spall A7-L5 

created under a pressure of p 0 =5.45 GPa. The 
digitized 3D view results from the traces 
obtained with a surface profilometer on the 
surface of the plastic replica. The projection on 
a plane perpendicular to the radial direction 
shows the height iso-contours in the spall. 
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Table 5.13 Spalls Geometry 


Samp. 

ID 

Track 

ID 

Hard. 

Imp. 

Ball 

Lot 

W^ax 

(mm) 

^max 

Wth 

kmax 

(mm) 

a 

o 

( ) 

Cut 

£ 

( ) 

D 

(mm) 

A1 

LI 

N 

R 

0.965 

1.20 

0.72 

62 

C 

21 

0.096 

A3 

LI 

Y 

R 

0.824 

1.02 

0.78 

54 

C 

21 

0.130 

A3 

L2 

Y 

R 

0.929 

1.15 

0.86 

62 

C 

20 

0.149 

A7 

LI 

N 

R 

0.988 

1.23 

0.99 

NA 

C 

NA 

0.138 

A7 

L3 

N 

R 

0.812 

1.00 

0.34 

NA 

C 

20 

0.075 

A7 

LI 

N 

R 

1.032 

1.28 

1.04 

30 

C 

24 

0.142 

A7 

L10 

N 

R 

1.082 

1.34 

0.91 

61 

C 

13 

0.130 

A1 

L2 

N 

S 

NA 

NA 

NA 

NA 

A 

54/68 

0.115 

A1 

L3 

N 

R 

0.747 

0.93 

0.57 

55 

A 

21/45 

0.089 

A3 

L3 

Y 

R 

0.845 

1.05 

0.85 

50 

A 

44 

0.142 

A7 

L4 

N 

R 

0.853 

1.06 

0.95 

30 

A 

67/65 

0.089 









A 

53/60 

0.130 

A7 

L5 

N 

R 

0.952 

1.18 

1.12 

45 

A 

10 

0.075 

A7 

L8 

N 

R 

0.900 

1.12 

1.07 

51 

A 

67/65 

0.108 


w max “ Maximum width of the spall, along axial direction (mm). 

Lmax " Maximum length of the spall, in the circumferential 
direction (mm) . 

a - Angle at the apex of the spall (°). 

/ 3 - Angle(s) observed in the respective cut (°). 

D - Depth, maximum depth for the circumferential cuts, and 
local depth for the axial cuts. 

C - Circumferential cut. 

A - Axial cut. 

N. A. - Not available 

All the tests included in this table were run at a peak Hertzian 
pressure of 5.45 GPa. 





CIRCUMFERENTIAL CUT 


RADIAL VIEW 



wmm 


' '////// 

AXIAL CUT 
(B-B) 


Figure 5.11 


Schematic representation of a spall indicating the 
nomenclature used in Table 5.13, angles are 
measured with respect to the surface of the 
sample . 
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the tip of the spall. Rough balls 
this test. 
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cantilever flap. Figures 5.13.a-d show top views of several 
spalls, as well as the result of axial sectioning. The axial cuts 
reveal the presence of subsurface cracks propagating in the 
axial direction. 

5.3 Retained Austenite ... 

The cyclic straining of the material under the contact may 
induce the gradual transformation of the retained austenite 
present in the matrix 26 . Test Series 2, see Table 5.2, sought to 
find a connection between the change in life and a probable 
change in retained austenite content. However, any changes in 
life are not obvious from the Table. More testing may be required 
before any conclusions may be drawn. 

5.4 Nucleation vs. Growth 

Experimental measurements summarized in Figure 5.14 27 show 
that the life of ground, 440C steel rods employed by the 3- 
ball-rod tester is reduced by about lOx when the conventional, 
lapped balls are replaced by the standard rough balls recommended 
for use with the testing machine. The installation of a 120 £im- 
diameter, 3 0 jxm-deep dent (Rockwell C hardness indent) reduces 
the contact life with rough balls by an additional factor of 3x 
to lOx. 

Since the effect of surface finish is confined to depths 
comparable to the asperity height, the implications of these 
results is that the lives of ground 440C bars obtained with the 
rough balls (standard conditions) and with smooth balls are 
controlled by the near surface mode. This is consistent with the 
findings of Lorosch shown in Figure 5.15 for the 52100 grade. 
This figure shows that only the elements with the longest lives, 
i.e. N » 10 6 contacts at p Q = 5.4 GPa, are of the subsurface 
mode . 

The installation of artificial dents makes it possible to 
follow the progress of spall nucleation and growth. Figures 
5. 16. a through c show the evolution of a surface initiated 
defect, in the proximity of a hardness impression. These pictures 
were taken from sample A6, subjected to p 0 =5.4 GPa, after 4, 6 
and 10.3 hours (approximately 2.1, 3.1 and 5.3 million cycles 
respectively). Figure 5.17 also shows the effect of a small 
surface indent in the raceway of a sample subjected to a rolling 
contact pressure of p 0 /k=8.9. 


26 The retained austenite may be as high as 45%. 

27 This figure is obtained from the results presented in 
Tables 5 . 5 (a) -5 . 12 . 
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Figure 5. 13. a Radial view of spall A3-L3 at 50X, the axial cut 

is shown at 80X. This experiment was run under a 
peak pressure of p 0 =5.45 GPa, with rough balls. A 
micro pit can be seen near the tip of the spall. 
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peak pressure of p 0 =5.45 GPa, with rough balls. 
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peak pressure of p 0 =5.45 GPa, with rough balls. 



peak contact pressure, MPa 


5000 


4000 


3000 


ZOOO 



1X10 5 lXlO 6 lXlO 7 1X10 8 lxlO 9 

stress cycles 


Figure 5.14 Fatigue lives versus contact stress level for 

440-C RCF test samples. Note: these are the 
nucleation and growth lives for unindented 
samples. The lives for dented samples are 
approximately equal to the growth part because 
nucleation is very short. The keys for the figure 
are indicated in the graph. 
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Figure 5.15 Cyclic lives and shakedown pressures for near- 
surface (C,D) ana subsurface rolling contact 
failures (A,B) in steel 52100 after Lorosch 
(1987) . Results labeled (C) are for 100 /xm surface 
indents . 
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Figure 5. 16. a Examples of early stages of spall formation near 
the hardness indent. The photomicrograph was taken 
after 4 hours of run time (approx. 2.1 million 
stress cycles) . The boxed region shown in (a) is 
magnified further in (b) . 






179 



$ in in 
* W *H 
C 

X! P 


-P ^ 
tfi rc 


C 
* 
o 
x o 

W *H 

i— I 

C/1 f — I 
•H *H 

fH 

(Q 

• (— i 

VO • 

»h r> 

in • 
X 

Q) o 

p p 

a a 

tP CL. 
*H fU 
U* w 

c ^ 

■h e 


U3 ~ 
C 2 

8 * 
^ P 
Q} 

0) +J 

X M-} 


c 

£ 

o 

JC 

w 


c 

o 

H 

Cn 

0) 

u 


X5 

0 ) c 

X -H 

o 

X P 

cu 

x: 

0) -P 

x: ^ 

Eh p 

4 -< 

. TJ 

in £ 
OJ-H 
rH c 

O CP 
>1 rc 

o e 


XI 

x> 

rH 

m 


<D— 


original page is 
of POOR QUALiTY 




180 




ORIGINAL PAGE IS 
Op POOR QUALITY 


Figure 5.16.C Fully formed spall near the location shown in the 
pryious two figures , after 10.3 hours (approx. 
5.3 million stress cycles). The boxed region shown 
m (a) is magnified further in (b) . 
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The growth life, N G , is defined as the portion of the total 
life consumed by the cyclic growth of a 10 jxm-surface crack into 
a full fledged spall, as seen in Figure 5.17. The corresponding 
initiation life, Nj, is obtained by subtracting N G from the total 
life, N, obtained by testing samples without a dent. The studies 
have shown that the growth life is relatively insensitive to the 
surface conditions and becomes a smaller fraction of the total as 
the severity of surface defects is reduced. In the case of the 
ground 44 0C steel rods tested at p 0 = 5.4 GPa, N G /N = 0.4 for 
rough balls and N G /N = 0.02 for smooth balls. The conclusion is 
that the cyclic growth of the spall can dominate the contact life 
of the near-surface initiated cracks when the surface is rough 
or damaged; nucleation of the spall dominates when the surface is 
smooth and undamaged. 

Preliminary studies show that nucleated cracks fail to grow 
at a contact pressure of p c = 24 GPa - This result points to the 
existence of a threshold for growth. 

5.5 Surface Roughness 

DeMul (1987) and Hamrock (1990), among others, have pointed 
out the important role played by surface asperities, bumps, 
grooves and dents in the near-surface failure mode. These defects 
cause stationary pressure spikes, conditioned by "micro" EHDL 
effects and produce peak amounts of plasticity at depths 
corresponding roughly with the surface regularity depth or height 
(2 to 2 0 jum) below the surface. 


Surface roughness features of the specimens and balls are 
characterized by peaks and valleys caused during the surface 
preparation. Figure 5.18 shows the running surface on specimen 
A6-L2 : a) after 0.26 million stress cycles, and b) after 2.73 

million cycles, c) spall formation after 10.09 million stress 
cycles. Examination of the running surface on the specimens 
indicates that small surface roughness features are gradually 
smoothed out with increasing numbers of stress cycles. However, 
deeper and more prominent features persist even after 15 million 
cycles. Similar observations are valid for the case of the 
roughened balls used in the RCF tests. The prominent peaks and 
valleys persist for millions of cycles and represent possible 
crack initiation sites. There is evidence (Figure 5.19) that 
nearby located microspalls may link together, thus providing the 
initial site for crack initiation. 


5.6 Replication and Metallographic Studies of Spalls in 440C 

There is a need for a more detailed description of the 3- 
dimensional geometry of the spalls, as well as of the different 
crack propagation mechanisms leading to that geometry. There have 
been previous three dimensional analyses of cracks under rolling 
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Figure 5.17 


Effects of a small, 120 /xm dent (A) in the raceway 
of a bearing steel RCF sample subjected to rolling 
contact at p 0 /k=8.9. (a) dent [A] before test, (b) 
crack nucleus [B] is visible on surface after N = 
0.62x10^ contacts, (c) N = 1.3xl0 6 contacts, and 
(d) spall [C] forms after N = 1.5xl0 6 contacts. 
The rolling direction is also indicated, i.e. R.D. 
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Figure 5.18 


Backscattered electron image of the specimen 
running surface (A6-L2) , after (a) 0.26xl0 6 

cycles, (b) 2.73xl0 6 cycles, and (c) lower 

magnification secondary electron micrograph of the 
spall formed after 2.73x1016 cycles. The boxed 
area shown in (c) is the location chosen for the 
micrograph in (a) and (b) . The rolling direction 
is also indicated, i.e. R.D. 
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Figure 5.19 Micrograph of the running surface of the 52100 RCF 

test balls showing two nearby located microspalls 
linked together. 
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contact by Kaneta and Murakami (1985a, 1985b, 1985c, 1986, 1987) . 
Although these address the problem of the lubricant induced 
opening mode and crack face friction, they failed to consider the 
cyclic plasticity induced under rolling contact loading. 

The top view and circumferential cuts shown in Figures 
5.12.a-d, reveal several important features: a) the cracks seem 
to have initiated at micro-pits near the tip of the spall (this 
is provided in the cases of Figures 5.12.b and c by the hardness 
indentation) , b) the crack at the V-apex slightly undercuts the 
surface, thus indicating a small amount of crack propagation in 
the reverse direction (more evident in 5. 12. a), c) the main 
cracks propagate at a shallow angle from the surface up to a 
depth of approximately 60-80 /m (evidence suggests that the crack 
front is straight along the axial direction and diverges from the 
initiation point at an angle of 50° - 60°), d) the cracks then 
propagate parallel to the surface for 130-200 /j,m , after which 
they change direction to an angle of 45° from the surface, e) 
this process is apparently being repeated cyclicly as evidenced 
by the unbroken flap still attached at the tip of the spall 
(clearly seen in 5. 12. a). 

The micrographs shown in Figures 5.13. a-d present axial cuts 
of different spalls. These figures suggest that subsurface crack 
propagation may also take place along the axial direction, as 
evidenced by the unbroken flap still attached to the surface. The 
results in Table 5.13 indicate that the maximum spall width is 
larger than the Hertzian (elastic) contact width, however the 
spalls were always confined to the wear track. Although some 
degree of conformity ("sink in") is expected to occur between the 
rod and ball, due to plastic deformation of the rod, the 
experimental evidence suggests that the balls are oscillating 
along the axial direction. This out-of-plane oscillation, and its 
associated shear tractions, may provide the necessary driving 
force for the subsurface, axial crack propagation. 

This cyclic change in directions would indicate that once the 
crack has been surface initiated by the micro pressure spikes, 
its propagation is the result of more than one mechanism. Bower 
(1988) has indicated that the fluid entrapment mechanism would 
play an important role in crack propagation. However this 
mechanism predicts that the crack would turn down, away from the 
surface, which is not supported by our experiments. 

The complicated non-proportional cycle of modes I and II 
stress intensities at the tip of the crack is more likely to be 
dictated by a combination of the entrapment mechanism and the 
resolved cyclic in-plane shear stress. Additionally, the cyclic 
bending to which the unbroken flap is subjected may favor the 
formation of cracks connecting the main crack to the surface. 

5.7 Rolling Contact of Hardened Aluminum 
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Rolling tests and finite element calculations using 
properties of 7075-T6 aluminum have been carried out as part of 
the present project. Aluminum was chosen primarily because the 
damage to the microstructure resulting in sub-surface crack 
formation is easily obtained within a short testing time, t = 0.5 
- 2.0 hours. In the case of 440C steel, much longer hours of 
operation are needed to cause crack initiation. Moreover, post- 
failure electron microscope studies of the spalls would indicate 
that they are primarily surface initiated cracks. 

Additionally, the cyclic stress-strain hysteresis response of 
7075-T6 aluminum is qualitatively similar to 440C bearing 
steels, as indicated in previous sections. Therefore, aluminum 
proves to be an ideal material to study sub-surface originated 
failure under rolling contact. 

The experiments and calculations performed on this material 
lead to the following conclusions: 

1) The depths at which cracks are observed to form in the RCF 
test samples are observed to correlate well with the maximum 
depth at which the maximum cyclic plasticity is calculated, see 
Figure 3.4.11. Subsurface crack initiation may play a more 
important role in this material. The accumulation of damage in 
the plastically strained region extending to a depth of 0.5w 
below the surface as shown in Figure 3.4.10, concurrent with a 
smoothing out of surface defects 28 , leads to the translating 
Hertzian pressure dominated failure. 

2) The region where the largest number of cracks is observed is 
the region where the residual stresses are not compressive and 
there is still some cyclic plasticity, see Figure 3.4.12. 
Compressive stresses inhibit cyclic crack growth by enhancing the 
frictional dissipation between the faces of the crack. 

5.8 Conclusions 

No connection was found between the amount of retained 
austenite and the life of the rolling element. More testing may 
be required before any conclusions may be drawn. 

The cyclic growth of the spall dominates the total life of 
the elements with surface initiated cracks, when the surface is 
rough or damaged. Nucleation is the controlling event in those 
samples with smooth and undamaged surfaces. 

Surface asperities, grooves, dents and bumps introduced 


28 This material shows strong evidence of lateral flow, i.e. 
ploughing, due to the rolling loading. 
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during manufacturing of the sample, play an important role in the 
near-surface failure mode. Prominent features are not smeared out 
even after several million cycles, thus providing sites for the 
highly localized, cyclic pressure spikes responsible for failure. 

The 3-dimensional features displayed by the spalls would 
indicate that the surface initiated cracks would propagate 
driven by several mechanisms. The individual contribution of the 
complicated, non-proportional cycle of Modes I, II and III, 
remains to be determined, and it certainly warrants further 
research in this area. 
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6. GENERAL DISCUSSION 


6.1 Finite Element Calculations of Contact Plasticity 

The contact plasticity has 2 components: (i) the transient 
part which is responsible for the radial displacements and 
residual stresses and (ii) the steady state or continuing 
plasticity which leads to damage accumulation and crack 
nucleation. There is evidence that the ELKP-model, which 
confines the transient to the first contact cycle, does not 
describe the transient adequately. The residual stresses for 2- 
dimensional, line contact derived from finite element calcula- 
tions with ELKP-properties understate the peak residual stresses 
by a factor of w 3 and place the peak at a relative depth of z/w 
= 1.3 instead of at z/w = 0.8, where it is observed (Hahn et 
al., 1987). Recently, the authors have treated the 2-dimensional 
contact with the more sophisticated, non-linear, Mroz-type, 2- 
surface representation devised by McDowell (McDowell, 1985 , 
Howell, 1991). In this case the transient extends over 15 and 
possibly more contacts, and while the residual stresses peak at 
the correct depth, the extrapolated values of the residual 
stresses are «3 times larger than the actual values. The steady 
state cyclic plasticity obtained with the non-linear model agrees 
with the results obtained with ELKP-behavior . It should also be 
noted that the 3-dimensional contact plastic zone observed in the 
aluminum samples in Section 3.4 is in good accord with the ELKP- 
calculat ions . These findings support the view that the 
descriptions of the continuing cyclic plasticity derived from the 
finite element calculations with ELKP-properties are reliable. 
However, the cyclic constitutive relations require more work to 
provide reliable descriptions of the transient and the residual 
stresses . 

While a single set of ELKP-properties has been used so far 
to characterized hardened steel, it must be noted that the 
kinematic yield strength, ct k , is sensitive to microstructure. An 
ongoing study by the authors shows that a K decreases by about 45% 
at room temperature as the amount of retained austenite increases 
from 16% to 46% . The loss of cyclic strength is accompanied by a 
more than 10-fold increase in the cyclic life and both of these 
changes occur even though the retained austenite is transformed 
by the cyclic deformation. While increases in the rolling 
contact life with retained austenite have been reported for 52100 
steel, these could not be demonstrated in the present work on 
440C steel . 

As demonstrated by the 3 dimensional finite element 
calculations, the addition of a sliding component with heating, 
results in the introduction of. tensile surface residual stresses. 
The peak value of cyclic plasticity takes place at approximately 
the same depth as for 3-dimensional pure rolling, i.e. z/w\ = 
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0. 3 29 , with approximately the same magnitude. However, the 
sliding component causes a non-zero cyclic plastic component at 
the surface. The combined effect of the tensile residual stresses 
and the cyclic plasticity may lead to an exhaustion of ductility, 

1. e. damage accumulation, thus giving support to the concept of 
surface originated spall failure (see Section 6.2). 

6.2 Spall Nucleation 

The present study demonstrates that the spall nucleation 
component dominates the rolling contact life when the contact 
surfaces are smooth, defect free and well lubricated. The work 
also shows, in agreement with previous studies, that modest 
changes in the surface roughness of the counterface alter the 
contact life by wlO-fold. The relatively large indents affect a 
«60-fold reduction in the contact life. These findings are 
consistent with recent analyses of the effects of surface 
irregulaties , which show that submicron-size grooves or 
asperities can produce large stationary pressure spikes near the 
defect (Goglia et al., 1984; de Mul et al., 1987; Hamrock, 1990). 
Rolling-plus-sliding, with attending frictional tractions, also 
produces more plasticity closer to or at the surface. As 
illustrated in Section 3.5, even larger amounts of near-surface 
plasticity that compete with the peak subsurface values are 
generated when the sliding is accompanied by intense frictional 
heating. Either individually or in combination, surface irregu- 
larities, sliding and frictional heating can shift the region of 
peak plasticity and damage from the subsurface to the near- 
surface, thereby altering the mode of spall failure. Work must 
be done to treat the plasticity arising from surface 
irregular ites so that this major contribution can be factored 
into the analysis of nucleation. 

The large reduction in the contact life produced by small 
indents may provide a means for accelerating contact testing. 
Service life is frequently limited by indents from debris or 
tramp abrasive particles and natural asperities. Consequently, 
rolling contact tests might prove more meaningful, and be of 
shorter duration, if the test samples were fitted at the start 
with well designed and pedigreed surface defects. This concept 
is exploited in the 3-ball-rod test which is designed to be 
conducted with artificially roughened balls. To provide a 
further test of this idea, the correlation of the nucleation 
lives obtained with and without defects should be examined. The 
contribution of inclusion particles may be a complicating factor 
here. This is because the cyclically deformed volume produced 
by a small, isolated, surface irregularity is so small that the 


29 Where w^ is the semi-major half-contact width; this is 
equivalent to z/w 2 = 0.53, in terms of the semi-minor half-contact 
width. 
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likelihood of an interaction with an inclusion particle is 
greatly reduced. 

In order to perform quantitative treatments of the spall 
nucleation life, the accumulated damage produced by the 
continuing cyclic plasticity and the conditions for crack 
nucleation must be treated. The studies performed on aluminum 
samples illustrate the complicated nature of the subsurface 
nucleation process. The results indicate that neither the rates 
of nucleation nor the rate of crack growth correlate with the 
amounts of cyclic plasticity. The findings suggest that residual 
stresses may be playing an important role. The reasons for the 
«10 4 -fold shorter crack nucleation lives at the same relative 
contact pressure displayed by the hardened aluminum compared with 
steel are not clear. It is possible that the higher homologous 
temperature of the aluminum tests and resulting greater atom 
mobility may facilitate microstructural changes leading to crack 
nucleation. Since microstructural changes also play a role in 
the nucleation of failure in steel (Swahn et al., 1976), the 
stability of specific microstructural elements may affect 
nucleation . 

6.3 Growth and Spalling 

The present study shows that the spall growth life is in the 
range N = 10 6 to 10 7 contacts, and is relatively insensitive to 
the condition of the surface. Similar results have been reported 
by Shao et al. (1987) for carburized steel. The implication is 
that the spall growth component is a large fraction of the total 
life when nucleation is "easy" and the total life is short, e.g. 
N < 10 7 contacts, and a tiny fraction of the total life when nu- 
cleation is difficult and the total life is long, e.g., N « 10 10 
contacts. The tests also offer evidence of a threshold for crack 
growth, i.e., crack nuclei about 10 /xm-long on the surface failed 
to grow at a contact pressure p 0 = 2.4 GPa. 

The fracture mechanics analyses examined in Section 4 reveal 
that the Hertzian stresses, centrifugal stresses, thermal stres- 
ses and the pressure of fluid in the crack cavity all contribute 
to a complex supperposition of Mode I, Mode II and Mode III 
crack driving forces. However, the largest contribution appears 
to be the Mode I driving force resulting from the fluid pressure. 
Further, it appears that 2-dimensional models of the 3- 
dimensional crack overstate the driving force. 

To critically test these findings, fracture mechanics calcu- 
lations of the spall growth life were carried out for the condi- 
tions in the 3 -ball-rod testing machine, and compared with the 
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measurements of the 440C growth life 30 . The results are 
presented in Figure 6.2. They show that the fracture 
calculations predict both the spall growth lives, at least within 
an order of magnitude. They also predict the existence of a 
threshold at roughly the observed contact pressure. Considerable 
uncertainty is involved in the calculations since the da/dN 
characteristics of the 440C steel were not known and the details 
of the dimensional growth were not considered. Further, the Mode 
II driving force was neglected. Consistent with this, the 
calculated lives are longer than the measured ones. 

The metallographic observations of a step-wise growth 
suggest that the spall advances into the surface with the lateral 
sides of the spall attached. At a later stage the latteral sides 
rupture converting the spall into a loose flap. After 
propagating parallel to the surface for a short distance, the 
spall branches away from the surface. Simultaneously, a crack 
seems to initiate from the surface attending to the cyclic 
tensile stresses due to the bending of the loose flap (Figure 
5.12a). This sets the stage for the formation of the first spall 
fragment. With the lateral portions of the remaining branch 
still intact, this branch resumes growth into the material until 
the lateral connections again break to produce a flap. One 
important difference is that large, rising and falling pressure 
spikes are expected to form around the cavity left by the first 
spall fragment. These pressure spikes will contribute to the 
crack driving force and accelerate the pace of spall growth. 

6.4 Conclusions 

(i) The findings reported here, taken together with other 
recent work, support the view that the descriptions of 
continuing cyclic plasticity derived from finite element 
calculations with ELKP-properties are reliable. However, the 
treatment of the cyclic constitutive relations must be improved 
to obtain reliable descriptions of the transient and the residual 
stresses . 


30 The total number of stress cycles was calculated by 
introducing the Mode I crack driving force derived in Section 4.4 
(Equation 4.1) into a Paris' type equation, i.e. da/dN = A (aK) n . 
The constants A and n fit the stage II portion of the results 
obtained by Bamberger et al. (1982), for different bearing steels 
(refer to Figure 6.1 and Table 6.1). Numerical integration was 
carried out assuming an initial crack nucleus of, a Q = 10 jum, and 
a final crack length of, a f = 240 ^ m. Small departures from the 
linear, stage II, portion of the resistance curve lead to a large 
(3-4 orders of magnitude) increase in total number of stress 
cycles. This can be related to the pressence of an effective 
threshold to crack propagation. 
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(ii) Either individually or in combination, surface 
irregularities, rolling-plus sliding and frictional heating can 
shift the region of peak plasticity and damage from the 
subsurface to the near-surface, thereby promoting near-surface 
spall failure. 

(iii) The spall nucleation component dominates the rolling 
contact life of 440C steel when the contact surfaces are smooth, 
defect free and well lubricated. Small surface irregularities 
reduce the nucleation life drastically. 

(iv) The possibility of accelerating rolling contact testing 
by fitting samples at the start with well designed and pedigreed 
surface defects deserves further study. 

(v) Studies using hardened aluminum as a model material 
suggest that residual stresses, atom mobility and the stability 
of microstructural elements affect the nucleation life. 

(vi) The spall growth component is a large fraction of the 
contact life when nucleation is "easy" and the total life is 
short, e.g. N < 10 7 contacts, and it is a tiny fraction of the 
total life when nucleation is difficult and the total life is 
long, e.g. N « 10 10 contacts. 

(vii) The experiments provide evidence of a threshold for 
crack growth. 

(viii) The measurements of the spall growth life are in 
reasonable accord with fracture mechanics predictions based on 
the Mode I component of the crack driving force generated by 3- 
dimensional cracks with fluid pressure in the crack cavity. The 
calculations predict both the growth lives within an order of 
magnitude and the existence of a threshold close to the observed 
contact pressure. 
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Table 6.1 Heat Treaments and Fracture Toughness 
Data for Bearing Steels, after 
Bamberger et al. (1982). 


Steel 

Heat 

Treatment 

Hardness 

Rc 

AK t h 

(MPa.ym) 

Kic 

(MPa. 7m) 

M50 

A 

63 

3.1 

17.7 

M50 

B 

61 

4.2 

23.0 

M5 0 

C 

43 

5.0 

48.0 

M50 

D 

44 

3.8 

51.0 


Heat Treatment 


A Austenitized at 1095 °C 

Tempered at 540 °C, 3 times. 

B Austenitized at 1095 °C 

Tempered at 540 °C, 5 times. 

C Austenitized at 1095 °C 

Tempered at 650 °C, 5 times. 

D Austenitized at 1095 °C 

Cooled at -85 'C 
Tempered at 540 °C, 3 times. 



Peak Pressure [MPa] 
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6.2 Fatigue Life (number of stress cycles to failure) 

versus the Contact Stress (peak Hertzian pressure) 
for bearing steel 440C. The figure indicates the 
experimentally obtained growth lives obtained for 
tests run with rough and smooth balls. The band 
represents the upper and lower limits for the 
total number of cycles calculated using the 
results from Section 4.4, and the data from Figure 
6.1 and Table 6.1. 
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7. CONCLUSIONS 

The following general conclusions can be drawn from the work 
performed for the present contract: 

The ELKP-parameters derived from the hysteresis loops of the 
440C steel are similar to those previously reported for hardened 
AISI 1070 and 52100 bearing steels. After correction for strain 
rate effects, cr^ = 1050 MPa and M = 188 GPa. 

The resistance of aluminum 7075-T6 under cyclic torsion is 
1/3 to 1/2 of that displayed by 440C steel. However the 
hysteresis loops reveal a kinematic behavior analogous to the 
440C. 


Hardened aluminum proves to be a reasonably good model 
material. Cyclic plasticity can be revealed metallographically , 
subsurface failure takes place after a moderate number of stress 
cycles, even at low contact pressures. These qualities may 
provide opportunities for: i) studying the effects of residual 
stresses, ii) verification of constitutive formulations, iii) 
testing fracture mechanics approaches to spall growth analyses. 

Three dimensional finite element calculations of repeated 
frictionless rolling contact have been carried out for ELKP 
properties of an AISI 44 0C bearing steel, at relative Hertzian 
pressures of Po/k^ « 4.0, 6.6 and 8.9, using the geometry of the 
three ball/rod contact fatigue testing machine. The subsurface 
stresses, plastic strains, residual stresses and plastic work 
done are evaluated. Equivalent stresses above the kinematic yield 
strength extend to the surface only for the highest load level. 
The peak equivalent plastic strain occurs at approximately 0.25WJ 
(w 1 /w 2 = 1.77) for all load cases; the values of these strains 
are 4.4xl0 -4 , 2.33xl0 -3 , and 4.0xl0 -3 , for the low, medium and 
high loads, respectively. 

Three dimensional finite element calculations of repeated 
frictionless rolling contact have been carried out for ELKP 
properties of the 7075-T6 aluminum alloy, at a relative Hertzian 
pressure of p Q /kk ~ 5.7, using the geometry of the three ball/rod 
contact fatigue testing machine. The subsurface stresses, plastic 
strains, and residual stresses are evaluated. The peak plastic 
strain amplitude, AeP ma x/ 2 = 1.6xl0“ 3 , occurs at a relative depth 
of z/w = 0.4 below the surface. Peak residual stresses obtained 
are a r x (max) = 2 6.4 MPa, and a r y(max) = 39.7 MPa, occuring at a 
relative depth of z/w =0.9. ; 

Two dimensional finite element calculations of repeated 
rolling-plus-sliding contact have been carried out for ELKP 
properties of an AISI 440C bearing steel, at relative Hertzian 
pressures of Po/k^ = 5.0. Plastic strains, residual stresses and 
temperature distributions are presented. Plastic strains are 


197 


shown to be considerable lower for ELKP properties than they for 
an elastic-perf ectly-plastic material (Kulkarni et al . , 1989). 
Steep temperature gradients close to the surface produce high 
compressive stresses; the ensuing cooling causes residual tensile 
stresses, a possible mechanism for thermo-cracking. Thermo- 
mechanical behavior is highly dependent on the material 
properties, which also may be temperature dependent. A more 
precise definition of these properties, especially the thermal 
film coefficient, is recommended. 

Efforts have been made to extend the Bower treatment of 2- 
dimensional, surface breaking cracks with fluid in the crack 
cavity to smaller relative crack lengths but the Bower algorithm 
proved unstable for a/w < 0.5. In the range 0.5 < a/w < 1.5, the 
values of the normalized driving force, K/p 0 ./a, decrease with 
the relative crack length, a/w, for lubricated contact, but 
increase for dry contact. 

A compilation of AK-driving force values of surface 
breaking cracks has been assembled. This reveals that fluid in 
the crack cavity, centrifugal stresses, and thermal stresses add 
significantly to the crack driving force, particularly for 
relatively short cracks, a/w < 0.1. Fluid in the crack cavity is 
especially potent, producing relatively large Mode I driving 
force values. The driving forces generated by 3-dimensional 
cracks are ”20% of the values calculated for 2-dimensional 
cracks . 

As a result of the modest driving force generated by small, 
3-dimensional cracks, surface cracks as large as a '10 /jl m would 
not be expected to grow in idealized smooth, well lubricated 
bearings with a Mode I, = 5 MPa subjected to a Hertzian 
Pressure, p 0 = 2.4 GPa. This is consistent with experimental 
results described in Section 5.4. 

No connection was found between the amount of retained 
austenite and the life of the rolling element. More testing may 
be required before any conclusions may be drawn. 

The cyclic growth of the spall dominates the total life of 
the elements with surface initiated cracks, when the surface is 
rough or damaged. Nucleation is the controlling event in those 
samples with smooth and undamaged surfaces. 

Surface asperities, grooves, dents and bumps introduced 
during manufacturing of the sample play an important role in the 
near-surface failure mode. Prominent features are not smeared out 
even after several million cycles, thus providing sites for the 
highly localized, cyclic pressure spikes responsible for failure. 

The 3-dimensional features displayed by the spalls would 
indicate that the surface initiated cracks would propagate 
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driven by several mechanisms. The individual contribution of the 
complicated, non-proportional cycle of Modes I, II and III, 
remains to be determined, and it certainly warrants further 
research in this area. 

The findings reported here taken together with other recent 
work support the view that the descriptions of continuing cyclic 
plasticity derived from finite element calculations with ELKP- 
properties are reliable. However, the treatment of the cyclic 
constitutive relations must be improved to obtain reliable 
descriptions of the transient and residual stresses. 

Either individually or in combination, surface 
irregularities, rolling-plus sliding and frictional heating can 
shift the region of peak plasticity and damage from the 
subsurface to the near-surface, thereby promoting near-surface 
spall failure. 

The spall nucleation component dominates the rolling contact 
life of 440C steel when the contact surfaces are smooth, defect 
free and well lubricated. Small surface irregularities reduce the 
nucleation life drastically. 

The possibility of accelerating rolling contact testing by 
fitting samples at the start with well designed and pedigreed 
surface defects deserves further study. 

Studies using hardened aluminum as a model material suggest 
that residual stresses, atom mobility and the stability of 
microstructural elements affect the nucleation life. 

The spall growth component is a large fraction of the 
contact life when nucleation is "easy" and the total life is 
short, e.g. N < 10 7 contacts, and it is a tiny fraction of the 
total life when nucleation is difficult and the total life is 
long, e.g. N ~ lo-*- 0 contacts. 

The experiments provide evidence of a threshold for crack 
growth . 

The measurements of the spall growth life are in reasonable 
accord with fracture mechanics predictions based on the Mode I 
component of the crack driving force generated by 3 -dimensional 
cracks with fluid pressure in the crack cavity. The calculations 
predict both the growth lives within an order of magnitude and 
the existence of a threshold close to the observed contact 
pressure . 
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Frictionless 3-D contact is simulated by translating a Hertzian pressure 
distribution across the surface of a semi-infinite body. This semi-infinite body is 
represented by a 3-D finite element mesh with appropriate displacement bound- 
ary conditions on the non-surface sides of the mesh. The boundary nodes are elas- 
tically displaced for each translating increment using the Boussinesq’s solution to 
the problem of a concentrated force acting on the surface of a semi-infinite body. 
For this purpose, the Hertzian pressure distribution is discretized into numerous 
concentrated forces. The boundary displacements due to each concentrated force 
are calculated and superposed for each increment. This appendix presents the 
detailed derivation of the Boussinesq’s solution for the case of the concentrated 
normal force, and then is extended for the case of the normal plus the tangential 
forces acting on the half-space. 

There are two ways of solving problems in elasticity. The first way is to 
find a set of stresses and strains which satisfy certain equations. This is the First 
Boundary Value Problem of Elasticity. The stresses must satisfy the equilibrium 
equations and the boundary conditions. The strains must satisfy the compatibility 
equations and the stresses must be related to the strains through appropriate 
stress-strain relations [See Appendix 2]. Alternately, as for Boussinesq’s solution, 
it is possible to reduce the above set of equations to three equations, the Navier 
equations, involving only the displacements, m. The problem is then known as the 
Second Boundary-Value Problem of Elasticity. The Navier equations in tensor 

notation are 


+ (X+m) v } { 


(i) 
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for simply connected domains bounded by a smooth surface, r where is the 
Laplacian operator, X and ^ or G are Lame’s constants, and v is the invariant 
defined as, v= u t t -. The solution of the Navier equations appropriate to the 
deformation of an elastic body by the concentrated force, K 0 , applied at some 
point, can be deduced from the particular integral due to Lord Kelvin, 


F(0 

u t ( x )= A f r [B— 



X 

3 




(2) 


where 


A= 

B= 


X-f fi 
X+3/z 

7 




and 


is the distance from the field point, X-, to the variable point, £. in r. It is sup- 
posed that the body forces F- are distributed over some subregion, of r, includ- 


ing the point, and vanish over the rest of the region. The resultant of the body 


forces acting on r 1 is 


F°=f F dr 

i J r t 


Let Fj increase in such a way that this integral has a finite limit, F^ 0 as ap- 
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proaches 0, to arrive at the notion of the concentrated force, F® acting at the 
point, Therefore, the displacements, u^x), produced at the point x i 7^ are 

F ° „ (x ■— £ )(x ( ■) 

X+ 3 /x t X+m t t j j_ po (3 

' Snfi{\+2n) r 87 t//(X+ 2^) r 3 J 


These expressions satisfy the homogeneous Navier equations at all points of the 
region except at the point of application of the force, x i = £;■ But if this point is 
deleted from the region by enclosing it in a sphere S of small radius a, the above 
solution in the remaining region corresponds to the deformation present in a 
body, r, with a cavity, S, subjected to the action of forces with the resultant, F°. 
For the case considered, let F-° act at the origin, and Fj Fg 0, and 

F 3 °=P, to get 


and 


u 


3 


=C{ 


X+3 fj 
(X+/*)r 


+ 



( 5 ) 


where C = — X+M , and r 2 —x x- . The stresses can be calculated by 

871>i(X+2 fi) 1 1 


r ij =xu k,k 5 lj + 
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The tractions produced by these stresses over the sphere S of radius r— a, are 
determined from 


where i /. = Xj/a. Integrating these tractions over the surface of the sphere, r=a, 
we obtain the components of the resultant force exerted on S by matter exterior 
to S. To solve the problem of deformation of the elastic half-space bounded by a 
plane subjected to the action of a concentrated force, it is necessary to account 
for the singularity at x ; = Boussinesq combined the earlier derived equations 
(11) and (12) with certain other singular solutions of Navier’s equations to get the 
dilationless (v=0) solution of Navier’s equations, so long as r 0, namely 

Dx , - 

a (6 

u =-7 r ' ' 

a r(r+ij) 

and 


D 


u 


3 


r 


( 7 ) 


where r 2 = x-X^ and D = constant. 

Now let us calculate the deformation of an elastic half space, x g > 0, 
due to a concentrated force, P, applied at the origin, acting in the positive direc- 
tion of the Xg-axis (see Fig. A.2.1) Since the point of application of the load is a 
singular point in the solution of Navier’s equations, delete it from x 3 > 0 by 
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describing a hemisphere of small radius 'sls. and consider the semi-infinite region 
bounded by the hemisphere and the x^-plane. Construct a solution such that 
the resultant of all external stresses acting on the hemisphere is P, and = 0, 
over the rest of the boundary. Adding equations (11) to (13), and (12) to (14) and 
applying appropriate boundary conditions, we get equations (15) and (16), 

„ Q = 1 , 2 ( 8 ) 

a 4717/ r 3 \+n r(i 3 +r) ; 

and 

u =—(— + ^ 1) (9) 

3 47r M l r 3 X +71 r 


Rewriting these equations using 

E 

2(l+i/) 


and 


Ev 


(l-f-^)(l— 2^) 

where E is the Young's modulus and is the Poisson’s ratio, we get 


wgr* _ (1 _ 2l/) 

a 2 nE v r 3 v 1 


r(z 3 +r) 


), a — 1,2 


( 10 ) 


and 
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(1+i/lP, 3 , . 1, 

“3=^f<7 + 2 ( 1 -'> ;> 


(ii) 


The same line of reasoning could be used to derive the appropriate 
boundary displacements for the case where the half-space is acted upon by a tan- 
gential force, T 1 . 


"x 2 ,E U 1 ' r (r+*f 2 ’ E V s 


{r+ Z y 


( 12 ) 


(i+i/)t,i y 2 „ w i y , , (i+^)-P r y^ 


r(r+z)‘ 


2ttE l f 3 




(r+z) J 


(13) 


(1 + v)T yz y I A ’ ~ 1 

^ = "i^"^ +(l ' 2l/) K^j ]+ ^ ^r + r 3^ 

The SUBROUTINE DISP in ABAQUS uses the last three equations to 
prescribe the elastic displacements on the non-surface boundary nodes for each 
translating increment. 


(i+*n 


.1 Z 


(14) 


3 It is related to the normal force P through the friction 
coefficient jx. 
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APPENDIX 2 

BOUNDARY CONDITIONS FOR 2-D FRICTIONAL CONTACT 


Frictional 2-D rolling contact is simulated by translating a Hertzian 
thermo-mechanical load across the surface of a semi-infinite body. The thermo- 
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mechanical load consists of three parts - 

X 2 j / 2 

1. normal mechanical pressure, p(x)=p o { 1 2 '’ 

W 

2. tangential surface traction, due to friction, f(x) = pp{x) and 

3. thermal load, due to the friction, q(x) = v t(x) 

The semi-infinite body is represented by a 2-D plane strain finite ele- 
ment mesh with appropriate boundaries. Describing the boundary conditions at 

the non-surface boundary nodes, consists of two problems, 

1. Prescribing the boundary temperatures due to the translating thermal 
load and 

2. Prescribing the boundary displacements due to the translating 
thermo-mechanical load. 

Boundary Temperatures 

Physical visualization of the problem, in conjunction with the work of 
Carslaw and Jaeger (1959), implies that the temperature gradients, which give 
rise to thermal stresses and distortions, are large only in the vicinity of the con- 
tact, extending to depth of a maximum of 0.5w. Hence, the temperatures at the 
bottom of the mesh could be fixed at the ambient temperature. In this part of the 
appendix, mathematical equations have been developed to reach the above con- 
clusion, which extend the work of Carslaw and Jaeger, to predict temperatures in 
a half space due to a translating non-uniform heat source. Carslaw and Jaeger 
derived the solution, for the temperature distribution in an infinite solid, due to 
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’Q’ units of heat/unit length liberated along the z-axis (See Fig. A.2.1), at the 
point (x’,0,y’). 

The temperature, T = t - * c , at a point (x,0,y) at time V is 


where 8 = ambient temperature, K = thermal conductivity, « — thermal dif- 
fusivity = K/pC, p = density and C = specific heat capacity . The result 
needed, for a half-space, can be obtained by appropriately choosing the heat 
source, Q, and integrating the above equation. The variation of the non-uniform 
translating heat source, distributed over a width 2w, is given by 


V(1 - l < 2 

W 


(15) 


and hence the limits of integration are from -w to w. Also, for temperatures at 
the current state, we have the heat source that has been moving for an infinite 
time at a velocity V across the half-space. Hence, the limits for the integration 
with respect to time are from 0 to oo, with the present time being oo. At an ear- 
lier time ’t’, the center of the heat source was located at -Vt. Now, the heat 
source ’Q’ is 2qdx’dt heat units/unit length (to account for the semi-infinite solid) 
through the point (x’-Vt,0,0). The temperature T’ is now given by 

adx'dt r ( x - x' + Vt ) 2 + y 2 , 

’ — 


T' 


2 nKt 


( 16 ) 
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The required result is obtained by integrating the above equation, using equation 
(20) with respect to ’x’ over the limits -w to w and with respect to time t from 0 
to oo. Hence, 


dt' 


T=I f t '=2 7k f-J 1 - J 1/2ix ' IT 


( x - £ + Vt ) 2 + y 2 
4ac t 


Evaluating the time integral, the above equation becomes 


Hf /-» K °<fj < 1 - > 2 + » 2 l 1/2 > - 4 ' ,W < 17 > 


where K (x) is the modified Bessel function of the second kind of order zero. To 
enable Gaussian integration the limits of the integral need to be changed to -1 to 
1. 

q —v [x — wu) 

T=- ^ /Lj e u K o {—[{x-wu) 2 W) 1/2 }{1 - u 2 ) l ! 2 wdu (18) 

For the calculations, v = 0.9 m/s, k — 1.4E-5 m^/s and the smallest value of the 
argument of K q comes out to be 192.0 and its value tends to zero and 6 tends to 
6 . Hence, the temperatures on the bottom of the mesh are fixed to ambient 
values and the effect of these temperatures in calculating the boundary displace- 


ments is neglected. 



Boundary Displacements 


In Appendix B, two methods of solving problems in elasticity were men- 
tioned and the second method, namely, the Second Boundary Value Problem of 
Elasticity, was used to derive the solution to Boussinesq’s problem. Here, for the 
2-D plane strain problem, the first method, namely, the First Boundary Value 
Problem of Elasticity will be used, in polar coordinates. The stresses must satisfy 
the equilibrium equations as well as the boundary conditions, the strains must 
satisfy the compatibility equations and the stresses must be related to the strains, 
by Hooke’s law. A suitable Airy Stress Function can be found to satisfy the equi- 
librium equations, the boundary conditons and the compatibility equations. From 
the Airy Stress Function, the stresses are calculated, then strains are found using 
Hooke’s law and the displacements are obtained from the strain-displacement 
relations and the boundary conditions. The derivations of the displacement equa- 
tions due to a concentrated vertical and horizontal line load indenting a half- 
space are presented in this section. 

a] Vertical or Normal line load at the surface of a half-plane : (see Fig. 

A.2.2). 

The appropriate Airy Stress Function is 

P 

JJ = r$s in 6 

7T 

1) The stresses are given by 
1 dU 1 d 2 U 


2 P cos 0 
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Then, using Hooke’s law, the strains are given by 

3 - k F(1 + *) cos 6 

2Ge rr ~ °rr ~ 4 ( a rr + ) = “ 2ir r 


2Ge 

— <7 

3 — /c 

rr 

rr 

4 

2G ^ 

b 

II 

= 0 

2Ge 00 

— a 66 

3 — K 
4 


( a rr + °ee ) — 


P [ 3 — k) cos i 


*'* t 6t~~ u 99 4 v'rr 1 2ir r 

But the strains are related to the displacements by the strain-displacement equa- 
tions, given by 

du 

2 Ge =2G—^ ( 22 

rr dr 

1, d “« 1 3 “r / 23 

2G <r , = 2G 7 |r-( T ) + ; ^-] l 23 


2Ge M = 2G ; lir + 

Using equations (23) and (26) and integrating, we get 

P 

2 G u = (1 + k) cos 0ln r + &(&) 

r 2 ^ 

Along the same lines, using equations (25) and (28), we get 

p P 

2Gu g = ^{l+ k) sin 9\n r - F[e) + g{r) + — (3 - «)sin 6 

Now, using equations (24), (27), (29) and (30), we get 

p ,9 o(r) 

F[6) + F" (fl) + - (ft - l)sin 0 + r — (— )=0 


A = F(<?) + F"(0) + - (is - l)sin 6 


-a = 4A 


then 


(29) 


I! II I 
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i E 

! fcJ 



The solution of the ordinary differential equation given by Eqn. (32) is 

p 

F[a) = C/ cos 6 + C 2 'sin 9 + A + — (k - l)<sin 0 

and the solution of the ordinary differential equation given by Eqn. (33) is 
g(r) = A + B'r, 

where A, B’, and C 2 ’ are constants. The displacements are found to be 

2 G u = [(1+k)cos 0ln r— ( k-I)cos fi+(/c-l)fein $]— C^sin 6+C ^ cos 6 

r 2jt 1 ' 

and 

2 Gu c = — [(l+»c)sin 0ln r— ( k— l)9cos 9+(3— >c)sin 9]— C^'cos 9— C 2 sin 6+ Br 
9 2w 

E 

For plane strain, 2 G — k — 3—4^ and hence we have 

u = ^il^! 4 (i/-l)cos 9ln r+2(l-2i/)(cos 9— 9sin ^J+CjSin e+C 2 cos S (30) 

T 27 tE 

and 

U — ^ 1+t ^i4(l-i/)sm flln r— 2(1— 2i/)9cos 9+4^sin flj+CjCos 9-C 2 sin 6+Br (31) 
0 2t tE 1 v 1 

The constants C p C 2 and B are calculated using the following boundary 

conditions - a) Points on the vertical axis have no lateral displacement, i.e. u^ = 

0 for 9 — 0. From Eqn. (34) 

Cj + Br = 0 i.e. Cj = B = 0 

b) For 9 = 0 and r = d, there is no vertical displacement. 

i-e- (u r ) fi=0 r=d = 0 

From Eqn. (34), we have 

(«,),_ = ^2(l-2,)]+C 2 

Hence, 

C 2 = -^^[4(i.-l)ln r+2(l-2i/)] 


and we get 
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u _ gH%,-l)cos 0ln (r/d)-2{l-2v)tem *)] 
r 27 tE 


and 


Ur = 


[4(1 — Asin flln (r/d)— 2(1— 2i/)flcos e+2sin 8 ] 
2 vE 1 v 


With Poisson’s ratio, «/ = 0.3, Young’s modulus, E - 207E9 Pa and 

( 1+0 


2nE 


and 


=9.9E-13 


u = p -X 9.9E— 13 [ - 2.8 cos 6 In (r/d) - 0.8 8 sin *)] 
r V 


(32) 


(33) 


= Py.x9.9E— 13[ 2.8 sin Sin (r/d) — 0.8 Scos S+2sin fl] 
rhese are the equations used in the SUBROUTINE DISP to prescribe the bound- 
L ry displacements due to the vertical or normal force, P y , by transforming these 
equations from polar to cartesian coordinate system, namely, 


and 


u = t^sin 8 + u e cos 8 


u = u cos 8 — u fl sin 8 
y r 


b] Horizontal force at the surface of the half plane : (see Fig. A.2.3) For 

this case the appropriate Airy Stress Function is also 
P . 

U — rSsin 8, 

n 

but the definition of 0 is different and from the earlier discussion we get the same 
equations (34) and (35). The constants C,, C 2 and B are different and come out 
to be : a) For 8 = tt/ 2 and r = d, there is no vertical displacement, i.e. 

KW/2,r=d = °‘ We haV6 ’ 


b) 


Cj=ir(l-2i/) 

Also, for 8=ir / 2 and r = d, there is no lateral displacement i.e. 
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( u e)fi=?r/ 2 . r: = d 
We have, 

4(1— i/)ln d-j-Av— Cy+Bd—Q 

3u r 

c) And, for 8 = tt/ 2 and r = d, (—) 6=jr / 2r=d = °* So ’ 

— A(u— l)ln d— 4(1— 2v)— C*2 = 0 

Solving these simultaneously, 

C 2 =-4(i/-l)ln d- 4(1-2j/) 

and 

B=~[4(*-l)] 

The equations now are 

u = 511^(4( I/ _1)C08 Sin (r/d)— 2(1— 2i/)(cos S+Ssin s)+*(l-2i/)sin s] 
f 2 ftE 


^ F(l+^ ) sta filn f J j d y 2 {\-2u)8cos s+4(l-i/>in 8 

v 2ft E 

4 r M 

+ 7r(l— 2i/)cos e+— (^-ljj 

Once again, - ' 0 9.995E— 13 And, 

27Tii 

u =P H X 9.9£-13[-2.8cos Sin (r/d)-0.8(cos S+Ssin S)+1.256sin S] (34) 


and _ 

2 8r 

u =P* X9.9£-13[2.8sin Sin (r/d)-O.80cos S+2.8sin S+1.256cos S —] (35) 

These equations are used in the SUBROUTINE DISP to prescribe the 

boundary displacements due to the horizontal force, P H , by transforming these 

equations from polar to cartesian coordinate system, namely, 
u - u.sin 8 — u cos S 

X v > 


u = u sin 8 + u,cos 8 
y r o 


and 


APPENDIX 3 


STRESS-STRAIN CONTOURS FOR 3-D ROLLING CONTACT IN AL 
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Some of the stress-strain contours calculated from the three 
dimensional finite element analyses of rolling contact in 
aluminum alloys are presented here. Figure A. 3.1 through A. 3. 6 
show the six components of stress when the contact pressure 
distribution is in the center of the mesh. Figures A. 3. 7 through 
A. 3. 16 illustrate the residual stress components and the residual 
plastic strain components after one contact sequence. The 
miniature schematic drawing by the side of each figure represents 
the sectional planes of the mesh on which the corresponding 
contours are calculated. 
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Figure A. 3.1 Axial stress (a x ) contours when the load is in 
the center of the mesh. 
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Figure A. 3. 2 Circumferential stress (ay) contours when the load 

is in the center of the mesh. 
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Figure A. 3. 3 Radial stress (a z ) contours when the load is in 
the center of the mesh. 
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Shear stress (^xy) contours when the load is in 
the center of themesh. 


Figure A. 3. 4 
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Shear stress (<7yZ) contours when the load is in 
the center of the mesh. 


Figure A. 3 . 5 
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Figure A. 3. 6 Shear stress (o xz ) contours when the load is in 

the center of the mesh.-.- 
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Figure A. 3. 11 Residual shear stress (cr r X z) contours 














Figure A. 3. 14 Radial residual plastic strain (e pr z ) contours. 
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Figure A. 3. 15 Residual shear plastic strain (tf pr X y) contours 



Figure A. 3. 16 Residual shear plastic strain (3’P r X z) contours. 




